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Abstract: We discuss various ideas on the nonperturbative vacuum struc- 
ture in QCD. The stochastic vacuum model of Dosch and Simonov is pre- 
sented in some detail. We show how this model produces confinement. The 
model incorporates the idea of the QCD vacuum acting like a dual supercon- 
ductor due to an effective chromomagnetic monopole condensate. We turn 
then to high energy, small momentum transfer hadron-hadron scattering. A 
field-theoretic formalism to treat these reactions is developed, where the ba- 
sic quantities governing the scattering amplitudes are correlation functions of 
light-like Wegner- Wilson lines and loops. The evaluation of these correlation 
functions with the help of the Minkowskian version of the stochastic vacuum 
model is discussed. A further surprising manifestation of the nontrivial vac- 
uum structure in QCD may be the production of anomalous soft photons in 
hadron-hadron collisions. We interpret these photons as being due to "syn- 
chrotron radiation from the vacuum". A duality argument leads us from 
there to the expectation of anomalous pieces proportional to (Q 2 ) 1 ^ 6 in the 
electric form factors of the nucleons for small Q 2 . Finally we sketch the idea 
that in the Drell-Yan reaction, where a quark- ant iquark pair annihilates with 
the production of a lepton pair, a "chromo dynamic Sokolov-Ternov effect" 
may be at work. This leads to a spin correlation of the qq pair, observable 
through the angular distribution of the lepton pair. 



1 Based on lectures given at the workshop "Topics in Field Theory" , organized by the 
Graduiertenkolleg Erlangen-Regensburg 1993 in Banz, Germany, and at the "35. Interna- 
tionale Universitatswochen fur Kern- und Teilchenphysik 1996" in Schladming, Austria. 
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1 Introduction 



In these lectures I would like to review some ideas on the way nonperturbative 
QCD may manifest itself in high energy collisions. Thus we will be concerned 
with strong interactions where we claim to know the fundamental Lagrangian 
for a long time now [|TJ: 

£qcd(x) = --G a Xp {x)G x ^{x) + £ q(x)(i^D x - m q )q(x). (1.1) 

Here q(x) are the quark fields for the various quark flavours (q = u, d, s, c, b, t) 
with masses m q . We denote the gluon potentials by G x (x) (a = 1, 8) and 
the gluon field strengths by 

G a Xp (x) = d x G a p (x) - d p G a x (x) - gf abc G b x (x)G c p (x), (1.2) 

where g is the strong coupling constant and f a bc are the structure constants 
of SU(3). The covariant derivative of the quark fields is 

D x q(x) = (d x + igG a x ^)q(x), (1.3) 

with A a the Gell-Mann matrizes of the SU(3) group. The gluon potential 
and field strength matrizes are defined as 

G x (x): = G a x (x)^, 

Gx P (x): = G a Xp (x)^. (1.4) 

The Lagrangian ( |1.1| ) is invariant under SU(3) gauge transformations. Let 
x — > U(x) be an arbitrary matrix function, where for fixed x the U(x) are 
SU(3) matrices: 

U(x)U\x) = 1, 

detU(x) = 1. (1.5) 

With the transformation laws: 

q(x) — > U(x)q(x), (1.6) 
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G x (x) -> U{x)G x {x)U\x) - ~U{x)d x U\x), (1.7) 

we find 

G Ap (x) -> f/(a;)G Ap (x)t/ t (x), 

and invariance of £qcd: 

£qcd(2:) -> £ QCD (x). (1.8) 



If we want to derive results from the Lagrangian ([LI), we face problems, 



the most notable being that £qcd is expressed in terms of quark and gluon 
fields whose quanta have not been observed as free particles. In the real 
world we observe only hadrons, colourless objects, quark and gluons are 
permanently confined. Nevertheless it has been possible in some cases to 
derive first principle results which can be compared with experiment, starting 
from £qcd ( |1 ■ 1| ) ■ These are in essence the following: 

(1) Pure short- distance phenomena: Due to asymptotic freedom the 
QCD coupling constant becomes small in this regime and one can make 
reliable perturbative calculations. Examples of pure short distance processes 
are for instance the total cross section for electron-positron annihilation into 
hadrons and the total hadronic decay rate of the Z-boson. 

(2) Pure long-distance phenomena: Here one is in the nonperturbative 
regime of QCD and one has to use numerical methods to obtain first principle 
results from £qcd, or rather from the lattice version of £qcd introduced by 
Wilson H . Today, Monte Carlo simulations of lattice QCD are a big industry 
among theorists. Typical quantities one can calculate in this way are hadron 
masses and other low energy hadron properties. (For an up-to-date account 
of these methods cf. ||). 

There is a third regime of hadronic phenomena, hadron-hadron collisions, 
which are - apart from very low energy collisions - neither pure long-distance 
nor pure short- distance phenomena. Thus, none of the above-mentioned 
theoretical methods apply directly. Traditionally one classifies high-energy 
hadron-hadron collisions as "hard" and "soft" ones: 

(3) High energy hadron-hadron collisions: 

(a) hard collisions, 

(b) soft collisions. 

A typical hard reaction is the Drell-Yan process, e.g. 
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tT + N -> 7* + X 

£+r (1.9) 

where £ = e, /i. All energies and momentum transfers are assumed to be 
large. However, the masses of the ir~ and N in the initial state stay fixed 
and thus we are not dealing with a pure short distance phenomenon. 

In the reaction ( |1.9|) we claim to see directly the fundamental quanta of 
the theory, the partons, i.e. the quarks and gluons, in action (cf. Fig. 1). In 
the usual theoretical framework for hard reactions, the QCD improved parton 
model (cf. e.g. ||), one describes the reaction of the partons, in the Drell-Yan 
case the qq annihilation into a virtual photon, by perturbation theory. This 
should be reliable, since the parton process involves only high energies and 
high momentum transfers. All the long distance physics due to the bound 
state nature of the hadrons is then lumped into parton distribution func- 
tions of the participating hadrons. This is called the factorization hypothesis , 
which after early investigations of soft initial and final state interactions || 
was formulated and studied in low orders of QCD perturbation theory in J7]]. 
Subsequently, great theoretical effort has gone into proving factorization in 
the framework of QCD perturbation theory [p|1-||10|. The result seems to be 
that factorization is most probably correct there (cf. the discussion in ||1 1|| ) . 
However, it is legitimate to ask if factorization is respected also by nonper- 
turbative effects. To my knowledge this question was first asked in [12] -[14 



In Sect. 4 of these lectures I will come back to this question and will argue 
that there may be evidence for a breakdown of factorization in the Drell-Yan 
reaction due to QCD vacuum effects. 

Let us consider now soft high energy collisions. A typical reaction in this 
class is proton-proton elastic scattering: 

p+p^p+p (i-io) 

at cm. energies E cm = a/s ~ 5GeV say and small momentum transfers 
\t\ = \q\ ^ lGeV. Here we have two scales, one staying finite, one going to 



infinity: 



E, 



cm 



<lGeV. (1.11) 
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Thus, none of the above mentioned calculational methods is directly appli- 
cable. Indeed, most theoretical papers dealing with reactions in this class 
develop and apply models which are partly older than QCD, partly QCD 
"motivated". Let me list some models for hadron-hadron elastic scattering 
at high energies: 



geometric | fL5|l , 



eikonal [16 



additive quark model |L7 
Regge poles fi~8H , 



topological expansions and strings |T9 
valons [pOfl , 

leading log summations pi 



two-gluon exchange ||22| , 

the Donnachie-Landshoff model for the "soft Pomeron' : 



It would be a forbidding task to collect all references in this field. The 
references given above should thus only be considered as representative ones. 
In addition I would like to mention the inspiring general field theoretic con- 
siderations for high energy scattering and particle production by Heisenberg 



JA\ and the impressive work by Cheng and Wu on high energy behaviour in 



field theories in the framework of perturbative calculations [25 



I will now argue that the theoretical description of measurable quantities 
of soft high energy reactions like the total cross sections should involve in 
an essential way nonperturbative QCD. To see this, consider massless pure 
gluon theory where all "hadrons" are massive glueballs. Then we know from 
the renormalization group analysis that the glueball masses must behave as 



m glueball OC 



Me -c/g 2 (M) 



;i.i2) 



for M — > oo, i.e. for g(M) — > 0, due to asymptotic freedom. Here M is the 
renormalization scale, g(M) is the QCD coupling strength at this scale and 
c is a constant: 



g\M) 



12 



4tt 2 331n(M 2 /A 2 ) 
8tt 2 



for M — > oo, 



11 



A : QCD scale parameter. 



;i.i3) 
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Masses in massless Yang-Mills theory are a purely nonperturbative phe- 
nomenon, due to "dimensional transmutation". Scattering of glueball-had- 
rons in massless pure gluon theory should look very similar to scattering of 
hadrons in the real world, with finite total cross sections, amplitudes with 
analytic t dependence etc. At least, this would be my expectation. If the 
total cross section cr tot has a finite limit as s — > oo we must have from the 
same renormalization group arguments: 

lim a tot (s) oc M"V C/92(M) (1.14) 

s— >oo 

for g(M) — > 0. In this case, the total cross sections in pure gluon theory 
are also nonperturbative objects! It is easy to see that this conclusion is not 
changed if cr tot (s) has a logarithmic behaviour with s for s —* oo, e.g. 

o-tot(s) -> const x (logs) 2 . (1.15) 

I would then expect that also in full QCD total cross sections are nonper- 
turbative objects, at least as far as hadrons made out of light quarks are 
concerned. 

Some time ago P.V. Landshoff and myself started to think about a possi- 
ble connection between the nontrivial vacuum structure of QCD - a typical 
nonperturbative phenomenon - and soft high energy reactions |[26|| . In the 
following I will first review some common folklore on the QCD vacuum and 
discuss in more detail the so-called "stochastic vacuum model". I will then 
sketch possible consequences of these ideas for high energy collisions. 



2 The QCD Vacuum 



According to current theoretical prejudice the vacuum state in QCD has a 
very complicated structure p7f-|37|. It was first noted by Savvidy |27] that 
by introducing a constant chromomagnetic field 



nr) a B, 



(2.1) 



into the perturbative vacuum one can lower the vacuum-energy density e(B). 
Here n and r] a are constant unit vectors in ordinary and colour space. The 
result of his one-loop calculation was 



e{B) = -B 2 + 
v ; 2 32tt 2 



In 



B 



(2.2) 
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where g is the strong coupling constant, M is again the renormalization scale, 
and (3q is given by the lowest order term in the Callan-Symanzik /3-function: 



For 3 colours and / flavours: 



/9b = 11 - \f- (2-4) 



Thus, as long as we have asymptotic freedom, i.e. for / < 16, the energy 
density e(B) looks as indicated schematically in Fig. 2 and has its minimum 
for B = -B vac 7^ 0. Therefore, we should expect the QCD-vacuum to develop 
spontaneously a chromomagnetic field, the situation being similar to that 
in a ferromagnet below the Curie temperature where we have spontaneous 
magnetization. 

Of course, the vacuum state in QCD has to be relativistically invariant 
and cannot have a preferred direction in ordinary space and colour space. 
What has been considered |53] are states composed of domains with random 



orientation of the gluon- field strength (Fig. 3). This is analogous to Weiss 
domains in a ferromagnet. The vacuum state should then be a suitable 
linear superposition of states with various domains and orientation of the 
fields inside the domains. This implies that the orientation of the fields in 
the domains as well as the boundaries of the domains will fluctuate. 

A very detailed picture for the QCD vacuum along these lines has been 
developed in ref. fl33"| . I cannot refrain from comparing this modern picture of 
the QCD vacuum (Fig. 4a) with the "modern picture" of the ether developed 
by Maxwell more than 100 years ago [|3l| (Fig. 4b). The analogy is quite 
striking and suggests to me that with time passing on we may also be able to 
find simpler views on the QCD vacuum. Remember that Einstein made great 
progress by eliminating the ether from electrodynamics. In the following we 



will adopt the picture of the QCD vacuum as developed in refs. [27|-[j53|],[i5S 
and outlined above clS cL working hypothesis . 

Let me now come to the values for the field strengths E a and B a in the 
vacuum. These must also be determined by A, the QCD scale parameter, 
the only dimensional parameter in QCD if we disregard the quark masses. 
Therefore, we must have on dimensional grounds for the renormalization 
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group invariant quantity (gB) 2 



(gBf oc A 4 . (2.5) 

But we have much more detailed information on the values of these field 
strengths due to the work of Shifman, Vainshtein, and Zakharov (SVZ), who 
introduced the gluon condensate and first estimated its value using sum rules 



for charmonium states 30 



<0|^G^(x)G^»(x)|0> = <0||^(B a (x)B a (x)-E a (x)E a (x))|0> 

= G 2 = (2.4 ± 1.1) ■ l(T 2 GeV 4 

= (335 - 430MeV) 4 . (2.6) 



Here we quote numerical values as given in the review [RUJ . A simple analysis 
shows that this implies 

< (% 2 B a (x)B a (x)|0 >= - < (% 2 E a (x)E a (x)|0 >= n 2 G 2 ~ (700MeV) 4 . 

(2.7) 



To prove eq. ( |2.7|) we note that Lorentz- and parity-invariance require the 
vacuum expectation value of the uncontracted product of two gluon field 
strengths to be of the form 

< ^ 2 G%{x)G h pa {x)\0 >= (g, p g ua - g^W^ (2.8) 



where G 2 is the same constant as in ( |2.6|) . Taking appropriate contractions 
leads to (fOjD and 0. 

We find that < 0|B a (x) B a (x)|0 > is positive, < 0|E a (x) E a (x)|0 > 
negative! This can happen because we are really considering products of 
field operators, normal-ordered with respect to the perturbative vacuum. 
The interpretation of eq. ( |2.7|) is, therefore, that the B-field fluctuates with 
bigger amplitude, the E-field with smaller amplitude than in the perturbative 
vacuum state. 

What about the size a of the colour domains and the fluctuation times r 
of the colour fields? On dimensional grounds we must have 

a~r~A _1 . (2.9) 
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A detailed model for the QCD vacuum incorporating the gluon condensate 
idea and a fall-off of the correlation of two field strengths with distance was 
proposed in pOJ: the "stochastic vacuum model" (SVM). In the following we 



will discuss the basic assumptions of the model and then apply it to derive the 
area law for the Wegner- Wilson loop, i.e. confinement of static quarks. For 
the rest of this section we will work in Euclidean space-time. To accomplish 
the analytic continuation from Minkowski to Euclidean space-time we have 
to make the following replacements for x x and G x (cf. (1.4)): 

x° — > —1X4, 

X, 
G° -> 

G^-Q. (2.10) 

Here X = (X, X 4 ) denotes an Euclidean space-time point and Q a (a = 
1, ...,4) the Euclidean gluon potential. With Q2.10|) we get 



(x-y)^-(X-Y) = -X a Y a , 

-ig J dx'G^x) -> -ig J dX a G a (X), 

G lk - g jk , (2.11) 

where 1 < j, k < 3, 1 < a, (3 < 4 and 

G a p = d a Gp - dpG a + ig[Q w Gp] (2-12) 
is the Euclidean gluon field strength tensor. 

2.1 Connectors 

Consider classical gluon fields in Euclidean space-time. Let X, Y be two 
points there and Cx a curve from X to Y (Fig. 5). We define the connector, 



the non-abelian generalization of the "Schwinger string" [[|l|] of QED as 



V(Y,X;C X ) = P{expH 5 f dZ a G a (Z)]}. (2.13) 

J C x 
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where P means path ordering. The connector can be obtained as solution of 
a differential equation. Let 

Z(n) = X, Z(t 2 ) = Y, (2.14) 

be a parametrization of Cx- Consider the differential equation for a 3 x 3 
matrix function V(r): 

±. V{r) = ^ g ^lg aiZ (r))V(r), (2.15) 
with the boundary condition 

V(n) = 11. (2.16) 

The solution of ( |2.15D , ( |2.16|) gives for r = r 2 just the connector Q2.13Q . 
Under a gauge transformation 

g a {x) -> u{x)g a {x)u\x) - l ~u{x)d a u\x) (2.17) 

9 

where U(X) G SU(3), we have 

V(Y, X; C x ) -> U(Y)V(Y, X; C X )XJ-\X). (2.18) 

The connector can be used to "shift" various objects from one space-time 
point to another in a gauge- covariant way. We define for instance the field 
strength tensor shifted from X to Y along Cx as 

g al3 (Y,X;C x ) :=V(Y,X;Cx)Gap(X)V~ 1 (Y,X;C x ). (2.19) 

Under a gauge transformation Qa^iY, X; Cx) transforms like a field strength 
tensor at Y: 

G aP (Y,X;Cx) -> U(Y)g a0 (X,X;Cx)U-\Y). (2.20) 

Connectors can, of course, be defined for arbitrary SU(3) representations, 
not only for the fundamental one used in (|2.13| ). Let T a (a = 1, ...8) be the 
generators of SU(3) in some arbitrary unitary representation R where 

[T a ,T b ] = if ahc T c . (2.21) 
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We define the connector for this representation by: 

V R (Y,X;C X ) : = P exp[-ip / dZ a G^Z)T a ). (2.22) 
We list some basic properties of connectors: 

(i) For 2 adjoining curves Ci,C 2 (Fig. 6), the connectors are multiplied: 

V R (X 3 ,X 1 ;C 2 + C 1 )=V R (X 3 ,X 2 ;C 2 ).V R (X 2 ,X 1 ;C 1 ). (2.23) 

(ii) Let C be a curve from X to Y and C the same curve but oriented in 
inverse direction, running from Y to X (Fig. 6). Then 

V R (X, Y; C)V R (Y, X; C) = 11, (2.24) 

Vl(Y,X;C) = V R (X,Y;C). (2.25) 

The product of the connectors for the path and the reverse path is equal to the 
unit matrix. The reversal of the path is equivalent to hermitian conjugation. 
We leave the proof of ( p.23| )-( ^25| ) as an exercise for the reader. 

2.2 The non-abelian Stokes theorem 

In this subsection we will derive the non-abelian generalization of the Stokes 



theorem Consider a surface S in Euclidean space time with boundary 
C = dS. Let X be some point on C as indicated in Fig. 7 and consider the 
connector ( |2.22| ) from X back to X along C: 

V R (X,X;C)=P exp[-ig f dZ a G a a (Z)T a }. (2.26) 

The problem is to transform this line integral into a surface integral. 

We start by considering a point Z\ in S and a small plaquette formed 
by curves C\, ..,64 where one corner point is Z\ (Fig. 7). We choose a 
coordinate system on S in the neighbourhood of Z\. 

{u,v)^Z(u,v) (2.27) 

such that 

Z 1 = Z(0,0), 
Z 2 = Z(Au,0), 
Z z = Z(Au,Av), 

Z 4 = Z(0,Av) (2.28) 
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and the curves Ci, Cz(C2, Ck) correspond to lines v =const. (it=const.). The 
matrix representing the line integral around the small plaquette is 

U R (Au, Av) := VRiC^VniC^Vni^VRid). (2.29) 

We want to make a Taylor expansion of U R (Au, Av) in Au and Av. From 
( |2.24| ) we find immediately that 

U R (0,Av) = U R (Au,0) = 1L (2.30) 

Thus the lowest order term in the expansion after the zeroth order is pro- 
portional to Am • Av and we get easily: 

U R (Au, Av) = 11 - ig X -Aa^G a ^{Z x )T a 

+0(Au 2 Av, AuAv 2 ), (2.31) 

where 

Au a /3 = AuAv d % a > Z e\ (2.32) 
o(u,v) 

In the limit Am, At> — > A<r al3 becomes the surface element da a p of the 
plaquette. 

The next step is to choose an arbitrary point Y, the reference point, on 
the surface 5* and to draw a fan-type net on S as a spider would do (Fig. 
8). The system of curves of the net consists of the curve Cx running from 
X to Y, then Cz 1 from Y to Z\, then around a small plaquette at Z x , back 
to Y along Cz 1 and so on. The final curve is Cx from Y to X. Apart from 
the initial and final curves Cx and Cx we have a system of plaquettes with 
"handles" connecting them to Y. With the help of Q2.23D -( |2.25| ) we see that 



the connector along the whole net is equivalent to the original connector 



V R (X, X; C) = V R (X, Y; Cx) ■ product of connectors 

for the plaquettes with handles • V R (Y, X; Cx)- (2.33) 

Let us consider one plaquette with handle, for instance the one at Z n in Fig. 
8. For this contribution to fl2.33p we get 

V R {Y, Z n - C z jV^(plaquette at Z n )V R (Z n , Y; C Zn ) 



V R (Y,Z n ;C Zn 



11 - ig-Aa aP g a aP {Z n )T a + 



V R (Z n ,Y;C Zn 



11 - ig l -Aa a ^{Y, Z n] C Zn )T a + ■■■. (2.34) 
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Here we use ( |2.31|) and the shifted field strengths as defined in ( [2.19| ). We 
leave it as an excercise to the reader to show that ( p?.19|) implies 



V R (Y, Z n - C Zn )g^(Z n )T a V R (Z n , Y; C Zn 
= g a (Y,Z ni C Zn )T a 



(2.35) 



for arbitrary representation R of SU(3). 

Inserting ( |2.34|) in (|2.33|) and summing up the contribution of all pla- 
quettes with handles, where we have of course to respect the ordering, we 
get in the limit that the net is infinitesimally fine: 



V R (X,X;C) = V R (X,Y;C X ) 



P exp 



^J s dcr aP g a ap (Y,Z;C z )T a 



V R (Y,X;C X ). (2.36) 



Here P denotes the ordering on the whole surface as implied by the net. 
Usually one takes the trace in ( |2.36| ) which leads with (|2.24]) to 



Tr V R (X,X;C) = Tr Pexp 



i g ~j s &cj a pg a aP (Y,z-c z )T a 



(2.37) 



This is the desired non-abelian version of Stokes' theorem. We leave it to the 
reader as an exercise to show that for the abelian case ( |2.37| ) reduces to the 
conventional Stokes theorem. 



2.3 The cumulant expansion 

As a last mathematical tool for making calculations with the SVM we discuss 



the cumulant expansion ||43|| . Consider functions 

t -> B(t) 



(2.38) 



on the interval O < r < 1 where B{r) are quadratic matrices. We assume 
that an averaging procedure over products of the functions B(-) is defined: 

E{B{ n )lE{B{ Tl )B{T 2 )),.... 

We consider first the case that all averages E(-) are c numbers and that 

£7(1) = 1. (2.39) 
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Let us consider the expectation value of the r-ordered exponential: 



f(t) :=£(Pexp 



t f 1 drBM 
Jo 



(2.40) 



where t e C. The cumulant expansion asserts that In f(t) can be expanded 



as 



oo ±n r \ „i 

ln/(t) = ^-/ dn.../ dr^^n, 
n! Jo Jo 



(2.41) 



where the n-th cumulant K n (ri, . . . , r n ) is a symmetric function of its argu- 
ments. A frequently used notation for the K n is: 



K n (r 1 ,...,r n ) = ((B(r 1 )...B(r n ))). 



(2.42) 



To prove ([2.41 ) we proceed as follows. Expanding in powers of t on the 
r.h.s. of ( |2.40| ) we get 

oo j.n 
fit) = 1 + E 



n=l 



nl 



where 



B n = f 1 d n . . . f 1 dT n E(P(B(n) . . . B(r n ))). 
Jo Jo 

Now we expand ln/(£): 



oo ±n 



71=1 



nl 



(2.43) 
(2.44) 
(2.45) 



where the expansion coefficients K n are obtained from: 



fit) 

oo f n 



1 + 



n=l 



nl 



exp(ln/(t)), 

OO ill 

exp ' x 



E rr^ 



n=l 



(2.46) 



From this we obtain the /C n as solution of the following system of equations: 



B 1 = K u 

B 2 = /C 2 + /Ci, 
3 

: /C 3 + -(/C 2 /Ci + /C1/C2) + /C 1? 

i3 4 = /C 4 + 2(/C 3 /Ci + /C1/C3) + 3/C2 + 2(/C 2 /C^ + /Ci/C 2 /Ci + K\K,<i) + /C^, 

(2.47) 



^2 
B 3 
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Clearly, this system can be inverted and we get K, n as sum of monomials of 
the form 

B h -B i2 ...B ik: (2.48) 

where 

5> = n. (2.49) 

Using Q2.44 ), every monomial (|2.48| ) can be written as n-fold integral over 
Ti, . . . ,r n with < Tj < 1 (j = 1, . . . , n). Since the integration domain is 
symmetric under arbitrary permutations of Ti, . . . , r n we can symmetrize the 
integrand completely. In this way we get 

K, n = [ dri . . . f dr n K n (ri, . . . , r n ), (2.50) 
Jo Jo 

where K n is a totally symmetric function of its arguments. Inserting Q2.50| ) 
in ( ^.45| ) we get the cumulant expansion ( |2.41| ), q.e.d. Explicitly we find for 
the first few cumulants: 

K x {\) = E{B{1)), 
K 2 (l,2) = E(P(B(1)B(2))) 

- 1 -E(B(1))E(B(2)) - 1 -E(B(2))E(B(1)), 

K 3 (l,2,3) = E(P(B(1)B(2)B{3))) 

~[E(P(B(1)B(2)))E(B(3)) 

+E(B(1))E(P(B(2)B(3))) 
+cycl. perm. ] 

+~[E(B(1))E(B(2))E(B(3)) + perm. ], 

(2.51) 

Here we write as a shorthand notation Ki(l) = Ki(ti), B(l) = B(ti) etc. 

The cumulant expansion (|2.41| ) has the so-called "cluster" property: Let 
us assume that the expectation values of the P-ordered products factorize 

E(P(B(1) . . . Bin))) = — {E(B(1)) . . . E(B(n)) + perm.} (2.52) 

n\ 
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for all n > 2 and all 

\n - Tj\ > Tmi n (i ^ j) ■ (2.53) 

We can then show that the cumulants for n > 2 vanish: 

tf n (l,...,n)=0, (n>2) (2.54) 

if the Tj satisfy ( |2.53| ). 

To prove ( |2.54|) we show first that it is true for n = 2. Indeed, from ( |2.51| ) 
and ( |2.52p we get for \t\ — r 2 \ > r min 

K 2 (l,2) = ±{E(B(l))E(B(2))+perm.} 

-\e{B{1))E(B{2)) - 1 -E(B(2))E(B(1)) 
= 0. (2.55) 

Now we proceed by mathematical induction. Assume that ( [2.54 ) has been 
shown for all k with 2 < k < n — 1. We have from Q2.471 ): 

K n {l,...,n) = E(P(B(l)...B(n))) 

-—[KM) . . .K x {n) +perm.l 

+ [symmetrized products of cumulants 2ffc(l, . . . , k) 
with 1 < k < n - 1] (2.56) 

but at least one factor with k > 2. With ( |2.52j ) we get now immediately 
K n (l, . . . , n) = in the region defined by ( |2.53|) , q.e.d. 

Up to now we have assumed the expectation values E(-) to be c-numbers. 
In this case many of the formulae ( [2.47| ) - ( |2.56|) can be simplified by using 

E(B(1))E(B(2)) = E(B(2))E(B(1)), 

etc. (2.57) 

We have, on purpose, not used such commutativity relations, since we are 
now going to generalize the cumulant expansion to the case where 

E(B(1)),E(B(1)B(2)),... (2.58) 
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are themselves quadratic matrix valued expectation values with 



£(11) = 11. (2.59) 

Then we have, of course, in general, no more the commutativity relations 
(|2.57 ). But all formulae ( |2.40 ) - (|2.56 ) are written in such a way that they 
remain true also for the case of matrix valued expectation values E(-). 

2.4 The basic assumptions of the stochastic vacuum 
model 

The basic object of the SVM is the correlator of two field strengths shifted 
to a common reference point. Let X, X' be two points in Euclidean space- 
time, Y a reference point and C x , C X i curves from X to Y and X' to Y, 
respectively (Fig. 5). We consider the shifted field strengths as defined in 
(|2.19|) and the vacuum expectation value of their product in the sense of 
Euclidean QFT: 



9 



2 



'Air 2 



1 



g^(Y,X;C x )g b pa (Y,X';C x ,) > =: ^F^X, X', Y; C x , C x 



(2.60) 

Here ( |2.20p and colour conservation allow us to write the r.h.s. of ( |2.6C| ) 
proportional to 5 ab . It is easy to see that F^ vpa depends only on X, X' 
and the curve C x + C x > connecting them; i.e. the reference point Y can 
be freely shifted on the connecting curve. In the SVM one makes now the 
strong assumption that the correlator ( j2.60| ) even does not depend on the 
connecting curve at all: 

• Ass. 1: F pupa is independent of the reference point Y and of the curves 
C x and C X i. 

Translational, 0(4)- and parity invariance require then the correlator 
( |2.60|) to be of the following form: 



d 

(Z a 5 w - Zpd^) + j^- (Z p 5 ua - Z a 5 up ) (1 - k)Di(-Z 2 )Y 

2.61) 



Fpupa — Ffj,vpa(Z) — ~^~*2\ {^pp^ua — ^pa^up) K-D(~Z 2 

+\\Jk {Za6 " p ~ ZAa)+ ifZp 
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Here Z = X — X', G2 is the gluon condensate, D, D\ are invariant functions 
normalized to 

D(0) = D x (0) = 1 (2.62) 

and k is a parameter measuring the non-abelian character of the correlator. 

Indeed, if we consider an abelian theory we have to replace the gluon field 
strengths by abelian field strengths J r llUl which satisfy the homogenous 
Maxwell equation 



^llVpo&v-F oa (^0 0, 



(2.63) 



if there are no magnetic monopoles present. It is easy to see that this implies 
k = in ( |2.61| ). Thus, in an abelian theory, the D-term in ( |2.61| ) is absent 
without magnetic monopoles but would be non-zero if the vacuum contained 
a magnetic monopole condensate. In the non-abelian theory the D-term has 
no reason to vanish. In fact, we will see that it dominates over the Dx-term. 
The abelian analogy suggests an interpretation of the D-term as being due 
to an effective chromomagnetic monopole condensate in the QCD vacuum. 
Two further assumptions are made in the SVM: 

• Ass. 2: The correlation functions D(—Z 2 ) and Di(—Z 2 ) fall off rapidly 

for Z 2 — > 00. There exists a characteristic finite correlation length a, 

which we define as ^ 

a:= ( dZ D(-Z 2 ). (2.64) 
Jo 

A typical ansatz for the function D, incorporating assumption 2, is as 
follows: 



D(-Z 2 ) 
which leads to 



27 
64 ( 



d*Ke lKZ K 2 



<aJ 



(2.65) 



D(—Z 2 ) oc exp 



3n\Z\ 



for Z £ 



00. 



(2.66) 



The function Di is chosen such that 



(2.67) 
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which leads to 

D^-Z 2 ) = (Z 2 )' 2 ( Z dv2vD(-v). (2.68) 
Jo 

With (p.67 ) the contracted field strength correlator has the form (cf. Q2.61|) : 

1 



flUflU 



:G,D(-Z l 



(2.69) 



The ansatz (|2.65|) , (|2.67|) can be compared to a lattice gauge theory calcu- 
lation of the gluon field strength correlator ffl] in order to fix the parameters. 
One finds (cf. 



45j and Fig. 9): 



a 

K 

G 2 



0.35 fm, 
0.74, 

(496 MeV) 4 , 



(2.70) 



with an educated guess for the error of ~ 10%. Of course, from Fig. 9 we 
get only the product k ■ G%. These quantities are obtained separately by 
measuring on the lattice several components of the correlator ( |2.60D . The 
value for G 2 in (|2.70|) from the lattice calculations is somewhat larger than 
from phenomenology. This can be understood as follows. The lattice calcu- 
lations of [[HJ are for the pure gluon theory. In the real world light quarks 
are present. Their effect is estimated to reduce the value of G 2 substantially 

We note that the correlation length is smaller, albeit not much smaller 
than a typical radius R of a light hadron (cf. e.g. P7| , |f48f , P9| , p]|j ) : 



Still 



R ~ 0.7 - 1 fm. 



a 2 /R 2 w 0.2 - 0.3 



(2.71) 



(2.72) 



is a reasonably small number and this will be important for us in the follow- 
ing. 

We come now to the third assumption made in the SVM: 



Ass. 3: Factorization of higher point gluon field strength correlators. 
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In detail assumption 3 reads as follows (cf. |45| ): 

All expectation values of an odd number of products of shifted field 
strengths vanish: 

(<?(l)...<?(2n + l)) = for n = 0, 1, 2, .... (2.73) 
Here we set as shorthand 

£(i) = ^i ft (Y,^;CxJ. (2.74) 
For even number of shifted field strengths we set in the SVM: 

(S(l)...S(2n)> = £ (G(h)G(j 1 ))...(g(i n )g(j n )), 

all pairings 

(ii,ji)--(i n Jn) (2.75) 

where n — 1,2, ... . 

We note that (£7(1)) must vanish due to colour conservation since the 
QCD vacuum has no preferred direction in colour space. The vanishing of 
the other correlators of odd numbers of field strengths, postulated in (|2.73|) , 



as well as the factorization property ( p. 75 ) are strong dynamical assumptions. 



They mean that the vacuum fluctuations are assumed to be of the simplest 
type: a Gaussian random process. 

For some applications of the SVM, for instance the calculation of the 
Wegner- Wilson loop described below, assumption 3 is not crucial and can be 
relaxed. But for the applications of the SVM to high energy scattering (cf. 
sect. 3) assumption 3 is crucial. In any case we prefer to specify the model 
completely, thus giving it maximal predictive power. On the other hand, of 
course, the model can then more easily run into difficulties in comparison 
with experiments. 



2.5 The Wegner-Wilson loop in the stochastic 
vacuum model 

We have now specified the SVM completely and can proceed to show how 
this model produces confinement . We consider a static quark-antiquark pair 
at distance R from each other and ask for the potential V(R). To calculate 
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V(R) we start with a rectangular Wegner- Wilson loop in the X\ — X4 plane 
(Fig. 10): Let C be the loop and S the rectangle, C = dS. Then 



W(C) = ^(Tr P exp[-ig j^Z^^Z)]) (2.76) 

and 

V{R) = -Y^\nW{C). (2.77) 

To evaluate W(C) in the SVM we first transform the line integral of the 
potentials in ( |2.76j ) into a surface integral of field strengths, using the non- 
abelian version of Stokes theorem (cf. sect. 2.2). For this we choose a 
reference point Y in S. We get then 

W(C) = i(Tr P exp[-ig f d^dX^Y, X, C x )]} (2.78) 

o J S 

where are the field strengths parallel-transported from X to Y along 
a straight line Cx- The path-ordered way to integrate over S in a fan- 
type net (cf. Fig. 8) is indicated by P. To evaluate the expectation value 
of the path-ordered exponential in ( [2 .78] ), we will use the technique of the 
cumulant expansion and the assumptions 1-3 of the SVM. We make the 
following replacements in the formulae ( |2.38| ) ff. of sect. 2.3: 



r ->(X ls X 4 ), 
B(Td^Gu(y,X®,C x v), 

E(B( Tl )..B(T n )) - ±Tr{g u (Y,xW ] C xW )...g u {Y,xW-,C x(n) )), 
t -> -ig, 

fit) - W(C). (2.79) 

With this we can express W(C) as an exponential of cumulants as given in 
(|2.41|). From Ass. 3 of the SVM we get that all cumulants for odd numbers 



of gluon field strengths vanish. The lowest nontrivial cumulant is the second 
one, K 2 , and from Q2.51 ) we find 



K 2 {1,2) -> ^Ti(P(g u (YX^;C xW )g u (Y,X^;C x(2) )). (2.80) 
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If we cut off the cumulant expansion Q2.41 ) at n = 2 we get then for the 
Wegner- Wilson loop 



W{C) =exp{-^- f dX l( iX 4 f dX[dX' A 
J s J s 

(\ti P [Qu(Y,X,C x )g u (Y,X',C x/ )])} 




jf dX lC lX 4 J s dX;dX^F 1414 (X - X') I , (2.81) 

where in the last step we used ( |2.61| ), i.e. Ass. 1 of the SVM. Next the 
assumption 2 of short-range correlation for the field strengths enters in a 
crucial way. This implies that for large Wegner- Wilson loops the integration 
over X[, X' A keeping X\, X 4 fixed in (|2.81| ) gives essentially a factor a 2 . The 
remaining X 1; X 4 integration gives then the area of S = RT. Thus we arrive 
at an area law for the Wegner- Wilson loop for R,T ^> a: 

W{C) = e~ uR - T , (2.82) 

where the constant a is obtained as 

^ f00 dZ*D(-Z*) = 327rKG * a2 . (2.83) 



a 



36 Jo 



We leave the proof of fl2.83j ) as an exercise for the reader. Comparing ( [2 .77]) 
and ( [2.82 ) we find that the SVM produces a linearly rising potential 



V{R) = a-R for R > a, (2.84) 

where a is the string tension. 

The results (2.82)-(2.84) were derived in the framework of the SVM in 
[ PH , where the model was introduced, and they are interesting in quite a 
number of respects: 

Only for k ^ does one get an area law and thus confinement. The D\ 
term which is the only one present in the abelian theory produces no con- 
finement. In the SVM confinement is related to an effective chromomagnetic 
monopole condensate in the vacuum. 

The short range correlation for the field strengths produces a long range 
correlation for the potentials if the D term is present in ( |2.61| ), i.e. if k ^ 0. 
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The string tension a is obtained numerically with the input ( |2.70| ) as 

a = (420 MeV) 2 (2.85) 

This is very consistent with the phenomenological determination of a ~ 
(430 MeV) 2 from the charmonium spectrum ]5T] . 



The SVM has been applied in many other studies of low energy hadronic 
phenomena (cf. |39j for a review). It was, for instance, possible to calculate 
flux distributions around a static quark- ant iquark pair |52fl . The results 



compare well with lattice gauge theory calculations wherever the latter are 
available. 

Let us come back to the calculation of the Wegner- Wilson loop above. 
It is legitimate to ask about the contribution of higher cumulants to W(C). 
How do they modify Q2.81Q ? It turns out that higher cumulants may cause 
some problems, which we discuss in Appendix A together with a proposal 
for their remedy. 



3 Soft Hadronic Reactions 
3.1 General considerations 

In this section we will present a microscopic approach towards hadron-hadron 
diffr active scattering (cf. J2f|, |53||). Consider as an example elastic scattering 
of two hadrons hx,h 2 

h 1 + h 2 ^h 1 + h 2 (3.1) 

at high energies and small momentum transfer. We will look at reaction 
( |3.1|) from the point of view of an observer living in the "femto-universe" , i.e. 
we imagine having a microscope with resolution much better than 1 fm for 
observing what happens during the collision. Of course, we should choose an 
appropriate resolution for our microscope. If we choose the resolution much 
too good, we will see too many details of the internal structure of the hadrons 
which are irrelevant for the reaction considered and we will miss the essential 



features. The same is true if the resolution is too poor. In ]53| we used a 
series of simple arguments based on the uncertainty relation to estimate this 
appropriate resolution. 

Let t = be the nominal collision time of the hadrons in (|3.1| ) in the 
cm. system. This is the time when the hadrons hi 2 have maximal spatial 
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overlap. Let furthermore be to/ 2 the time when, in an inelastic collision, 
the first produced hadrons appear. We estimate to ~ 2 fm from the LUND 
model of particle production |51 |. Then the appropriate resolution, i.e. the 



cutoff in transverse parton momenta kx of the hadronic wave functions to be 
chosen for describing reaction (|3.1| ) in an economical way is 



k 2 T < y/s/(2t ) (3.2) 

where \fs is the cm. energy. Modes with higher k^ can be assumed to be 
integrated out. With this we could argue that over the time interval 

- -to < t < -t (3.3) 

the following should hold or better: could be assumed: 

(a) The parton state of the hadrons does not change qualitatively, i.e. 
parton annihilation and production processes can be neglected for this time. 

(b) Partons travel in essence on straight lightlike world lines. 

(c) The partons undergo "soft" elastic scattering. 

The strategy is now to study first soft parton-parton scattering in the 
femto-universe. There, the relevant interaction will turn out to be mediated 
by the gluonic vacuum fluctuations. We have argued at length in section 2 
that these have a highly nonperturbative character. In this way the non- 
perturbative QCD vacuum structure will enter the picture for high energy 
soft hadronic reactions. Once we have solved the problem of parton-parton 
scattering we have to fold the partonic S'-matrix with the hadronic wave 
functions of the appropriate resolution ( ^.2| ) to get the hadronic S'-matrix 
elements. 

We will now give an outline of the various steps in this program. 

3.2 The functional integral approach to 
parton-parton scattering 

Consider first quark-quark scattering: 

q(pi) + q(p 2 ) -> q(p 3 ) + qfa), (3.4) 

where we set 
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s = (pi + p 2 ) 2 

t = (pl- P3) 2 

u=(Pi-Pa) 2 - (3.5) 

Of course, free quarks do not exist in QCD, but let us close our eyes to this 
at the moment. Now we should calculate the scattering of the quarks over 
the finite time interval ( j3.3|) of length to ~ 2 fm. Let us assume that 2 fm is 
nearly infinitely long on the scale of the femto universe and use the standard 
reduction formula, due to Lehmann, Symanzik, and Zimmermann, to relate 
the 5"-matrix element for ( |3.4|) to an integral over the 4-point function of the 
quark fields. We use the following normalization for our quark states 

(q(p j ,s j ,A j ) I q(p k ,s k ,A k )) 

= 5 SjSk 5 AjAk (2n)^2pppl5 3 ( Pj - p k ) 
= 6{j,k), (3.6) 
where Sj, s k are the spin and A,-, A k the colour indices. With this we get 

(q(p 3 , s 3 , A 3 )q(p A , s 4 , At)|%(pi, Si, A 1 )q(p 2 , s 2 , A 2 )) 
= (3,4|5|1,2> 

= (3,4(1, 2) + V {(4|(z? - m' q ) ® {2,\{i$ - m' q ) 
(0|T(g(4)g(3)g(l)g(2))|0) 
{i$+m' q )\\)®(0+m' q )\2)y (3.7) 

Here is the quark wave function renormalization constant and m' q the 
renormalized quark mass. We use a shorthand notation 

Q{j) = ?( x i)' 

(j = l,..,4), (3.8) 

where u is the spinor in Dirac and colour space. Two repeated arguments 
j, k, ... imply a space-time integration, for instance 

q{l){i$ + m' q )\l) = J dxtqixJM + m' q )e- ipixl u sl)Al { Pl ). (3.9) 
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Thus in ( |3.7| ) we have four integrations over xi, ...,x&. 

We can represent the 4-point function of the quark fields as a functional 
integral: 

(0|T(g(4)g(3)g(l)g(2))|0) 

= Z- 1 J V(G, q, q) exp {i J dx£ QCD (x)}g(4)g(3)g(l)<?(2), (3.10) 

where Z is the partition function: 

Z = J V(G,q,q)exp{i J dxC QC D(x)}. (3.11) 

The QCD Lagrangian fll.lj ) is bilinear in the quark and antiquark fields. 
Thus - as is well known - the functional integration over q and q can be 
carried out immediately. After some standard manipulations we arrive at 
the following expression: 

(0iT(g(4)g(3)g(l)g-(2))|0) 

= ly V(G)ex P l [ -t J d^Tr(G Ap (x)G A '(x))} 
det[— i(i^ x D\ — m q + it)} 

{t^(4, 2; G)^S F (3, 1; G) — (3 <-> 4) j . (3.12) 

Here Sp(j, k] G) = Sp(xj, x^, G) is the unrenormalized quark propagator in 
the given gluon potential G\(x). We have 

(i-fDp - m q )S F (x, y- G) = -6{x - y). (3.13) 

Functional integrals as in ( |3.12| ) will occur frequently further on. Let F(G) 
be some functional of the gluon potentials. We will denote the functional 
integral over F(G) by brackets (F(G))c- 

(F(G)) G := || X?(G)exp j-* J d^Tr(G Ap G A ')} • 
det[-i(i-f X D x -m q + ie)]F(G). (3.14) 



Now we insert ( |3.12| ) in ( J3.7\ ) and get 

7j i> 



(3,41511,2) = (3,4|1,2) - ZT 2 (M^(G)MUG) - (3 <-> 4)) G , (3.15) 
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where 



M^G) = (k\(i0 -m' q )S F (i0 +m' q )\j), 

(A; = 3, 4; j = 1,2). (3.16) 

The term A4£ • -Mf 2 on the r.h.s. of (|3.15|) corresponds to the t-channel 
exchange diagrams, the second term, where the role of the quarks 3 and 4 
is interchanged, to the w-channel exchange diagrams (Fig. 11). The latter 
term should be unimportant for high energy, small \t\ scattering. Thus we 
neglect it in the following. For the scattering of different quark flavours it is 
absent anyway. We set, therefore: 

(3, 41511, 2) = (3, 4|1, 2) - Z^{MUG)MUG)) g . (3.17) 

We can interpret Jid^JG) as scattering amplitude for quark j going to k 
in the fixed gluon potential G\(x). To see this, let us define the wave function 

\^) = S F (0+m' q )\j) (3.18) 
which satisfies the Dirac equation with the gluon potential G\(x)\ 

(i 7 A D A -m g )|V£) = 0. (3-19) 
Furthermore we use the Lippmann-Schwinger equation for S F : 

S F = SP - SP (g0 - 5m q )S F (3.20) 

where SP is the free quark propagator for mass m' q and 5m q = m' q — m q is 
the quark mass shift. Inserting ( 3.20j) and (|3.18 ) in (|3 .16 ) gives after some 



simple algebra 



M F k AG) = ( Pk \(gG-8mM). (3.21) 



This represents M-kj m ^ ne f° rm a scattering amplitude should have: a com- 
plete incoming wave is folded with the potential and the free outgoing wave. 
However, there is a small problem. The wave function |V^-) defined in ( p,18j ) 
does not satisfy the boundary condition which we should have for using it 
in the scattering amplitude, i.e. it does not go to a free incoming wave for 
time t — > — oo. The wave function with this boundary condition is obtained 
by replacing the Feynman propagator Sf in (|3.18| ) by the retarded 
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We have shown in [ jj^] that in the high energy limit this replacement can in- 
deed be justified for the calculation of Jv[^-(G) if the gluon potential G\(x) 
contains only a limited range of frequencies: 



M F kj {G) 



MUG) = (p k \(g0 - 5m)\if>. 



Pi- 



where 



which satisfies 



ivy 



S r (i0+m')\j), 



(i'fD lt -m q )\if>. 

for 



0. 



-oo. 



(3.22) 

(3.23) 

(3.24) 
(3.25) 

We summarize the results of this subsection: At high energies and small 
\t\ the quark-quark scattering amplitude can be obtained by calculating first 
the scattering of quark 1 going to 3 and 2 going to 4 in the same fixed gluon 
potential G\(x). Let the corresponding scattering amplitudes be Ai^G) 
and M. r A2 {G) (cf. ( |3.22j )). Then the product of these two amplitudes is to be 
integrated over all gluon potentials with the measure given by the functional 
integral ( |3.14| ) and this gives the quark-quark scattering amplitude (|3.17|) . 
The point of our further strategy is to continue making suitable high energy 
approximations in the integrand of this functional integral which will be 
evaluated finally using the methods of the stochastic vacuum model. 

In the following it will be convenient to choose a coordinate system for 



the description of reaction Q3.4|) where the quarks 1,3 move with high velocity 



in essence in positive x direction, the quarks 2,4 in negative x direction. 
We define the light cone coordinates 



X± = X ± X 



(3.26) 



and in a similar way the ± components of any 4-vector. With this we have 
for the 4-momenta of our quarks: 



1 \Vr 



+ 



2p j+ 



2P. 



(3.27) 



'3 + 
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for j = 1,3 with pj + — > oo and 

/ i 

Pk = 



Pk- 



for k = 2, 4 with — > oo. 



PfcT 



(3.28) 



3.3 The Eikonal expansion 



The problem is now to solve the Dirac equation ( 3.24 ) for arbitrary external 
gluon potential G\{x). Of course, we cannot do this exactly. But we are 
only interested in the high energy, small \t\ limit. This suggests to use an 
eikonal type approach. This works indeed, but it is not as straightforward 
as one would think at first, since the Dirac equation is of first order in the 
derivatives, whereas the eikonal expansion is easy to make for for a second- 
order differential equation. What we did in |53| was to make an ansatz for 
the Dirac field ip r p .{x) in terms of a "potential" <f>j(x) as follows: 

%.(x) = (i 1 x D x + m q )cj )j (x). (3.29) 

A suitable boundary condition for <fij(x) which is compatible with ( p. 25 ) is 

1 1 + 7° 



[x 



+ m q 2 



-i Pj 



for t — > —oo. Inserting (|3.29| ) into the Dirac equation ( p. 24]) gives: 



{z 7 A D A - m q }{iYD p + mJ^Or) = 0. 



(3.30) 



(3.31) 



For the case of no gluon field, G\(x) = 0, the covariant derivatives D x de- 



generate to the ordinary ones, d\, and (|3.31| ) to the Klein-Gordon equation: 

(D + ml)(f> j (x) = 0. (3.32) 



Thus the problem for the Dirac field is in essence reduced to a scalar field-type 
problem which is handled more easily. 
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Now it is more or less straightforward to turn the theoretical crank and 
to obtain the eikonal approximations for <j)j(x) and ip r p . (x), respectively. Take 
j = 1 as an example. We make the ansatz 

X ( X ) = e-^Mx). (3.33) 

Inserting pi from ( |3.27| ) we see that on the r.h.s. of ( |3.33| ) the fast varying 
phase factor is 

ex P(-«2^ 1 + a; -)' 

since pi + — > oo. Assuming all the remaining factors to vary slowly with x 
we insert ( |3.33| ) in ( |3.31| ) and order the resulting terms in powers of l/pi+. 
The solution of (|3.31|) to leading order in l/p\+ is then easily obtained. The 
final formula for (x) reads: 



P Pl {x) = V-(x + ,x-,x T ).{l + 0[ — 



-ipix 



«(Pl), 



where 



V-(x+, X-,x.t) = P \ exp 



:9 



dx' + G-(x' + , x_, x T ) 



(3.34) 



(3.35) 



and P means path ordering. When coming in, the quark picks up a non-abeli- 
an phase factor , just the ordered integral of G along the path. Of course, V_ 
is a connector, as studied in section 2.1, but now in Minkowski space and for 
a straight light-like line running from — oo to x (Fig. 12). 

In a similar way we obtain for the quark with initial momentum p 2 , i.e. 
the one coming in from the right, for p 2 _ — > oo: 

1 



r P2 (x) = 

where 

V + (x + , x_,x r 
Recall that 



V + (x + ,x_,x r ) jl + I — 1 \e-^ x u(p 2 



exp 



~2 9 



dx'_G + (x + , x'_,y. T ) 



(3.36) 

(3.37) 
(3.38) 



G±{x) = (G 0a (x)±G 3a (x)) y 

are matrices in colour space. Thus path ordering in ( 3.35Q and (|3.37|) is 
essential. 

A solution for ty T {ip r p2 ) as a series expansion in powers of l/pi + (l/p 2 _) 
was obtained to all orders in [55[|. 
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3.4 The quark-quark scattering amplitude 

We can now insert our high energy approximations ( |3.34j ) , ( |3.36| ) for ij) r \x) 



in the expression for Mt- in ( |3.22|) . The resulting integrals are easily done 
and we get for Pi+,P3+ — > oo: 



M r 31 (G) 



dxe i(p3 - pi)x 
u{pz){g0{x) - 5m)V-(x)u{pi) 



:(Pa ~ Pi)+x- - i(p 3 - Pi)t ■ x T 



— > —J dx + dx_d 2 Xx exp 
d 

u(p3h+-^-V^( y x + , x-, x T )«(pi) 
VP3+P1+ ' 8*3,81 J dx_d 2 x T 

exp^ (Pa - Pi)+a;- - ^(Ps - Pi)t ■ x T ] 
[y_(oo,x„,x T )-ll] A3iAl . (3.39) 

In a similar way we obtain for p2-,f>4- — ► 00: 

^42( G ) ■ P2- ■ Ss 4 ,S2 

J dy + d 2 ?/ T exp[^(p4 ~P2)-y+ - ~ Y>2)t ■ yr] 
[F + (|/ + ,oo,y T )-]lU A . (3.40) 

Here we have written out the spin and colour indices of the in- and outgoing 
quarks (cf. Q3.6Q , ( p.7|) ). We have used furthermore: 



UsAPih^Us^) 



VP3+P1+ ■ ^3,«i n + fOT Ply 



*2 



n_ 



for p 2 ,4- 



00, 

00, (3.41) 



where 



n± 



(3.42) 



( 1 \ 





V ±1 / 

Finally, the x + (x_) integration for -M^-M^) could be done with the help 
of fl2.15|) or rather the analogous equation for connectors in Minkowski space 
time. 
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Now we can insert everything in our expression ( |3.17 ) for the S'-matrix 
element. This gives: 

(3,4|5|1,2) = (3,4|1,2) + 

VP3+P1+P4-P2- • S S3)S1 S SA>S2 Z^ 2 J dx_d 2 x T J &y+& 2 y T 
exp [-(p 3 - Pi)+x- - i(P3 - Pi)t • x T ] 



ex P [g (P4 - P2)-y+ - Kpa - P2)t • yr] 
([y_(oo,x_,x r ) - nj^JV+^+jOOjyr) - U]a 4 ,a 2 ) g 

From translational invariance of the functional integral we have: 



(3.43) 



([V-{oo,X-,x T ) - n] As>Al [V + (y + ,oo,y T ) - U]a 4 ,a 2 )g 

= ([F_(oo,0,x r -y T ) - U] A3 , Al [K + (0,oo,0) - 11]a 4 ,a 2 > G . (3.44) 



Inserting ( |3.44 ) in Q3.43 ), we can pull out the 5-function for the overall energy- 
momentum conservation and we get finally: 

(3,4|S|1,2> = (3,4|l,2) + +2(27r) 4 5(p 3 +P4-pi-p 2 )(3,4|T|l,2), 

(3,4|T|1,2) = • Ss 3 , Sl S S4 ,s 2 (-Z-; 2 ) J d 2 z T e ic ™ 

([y_(oo, 0, z T ) - ii] A3iAl [y + (o, 00, 0) - nu^ a ) G . (3.45) 

Here g is the momentum transfer: 



9 = Pi-P3=Pi- Pa, 
g 2 = f. 

In the high energy limit g is purely transverse 

/ 



(3.46) 



q 




q 



-q|. 



(3.47) 



Using different techniques the type of formula ( p.45 ) was also obtained in 
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In ( |3.45| ) we still have to calculate the wave function renormalization 
constant Z^. This can be done by considering a suitable matrix element of 
the baryon number current 

-q(x)^q(x) 

which is conserved and, therefore, needs no renormalization. The result is 
(cf. 



-(TrV_(oo,0,0)) G . 



(3.48) 



Let us summarize the results obtained so far: 

The quark-quark scattering amplitude ( |3.45|) is diagonal in the spin in- 
dices. Thus we get helicity conservation in high energy quark-quark scatter- 
ing. Using ( [3.41|) we can write the spin factor in (|3.45| ) as 



2VP3+Pl+P4-P2-5 S3 ,si^4,82 - -S3 



{Psh fl U sl {pi)u SA (p4)juU S2 {p 2 



(3.49) 



for pi 5 3+ — ► oo and ^2,4- — * oo. This 7 M £g> 7 M structure was postulated for 
high energy quark-quark scattering in the Donnachie-Landshoff model for the 
"Pomeron" coupling ||23|| . However, a study of quark- ant iquark scattering (cf. 
[ [53H and below) reveals that (|3.45f) does not allow an interpretation in terms 
of an effective Lorentz vector exchange between the quarks. The amplitude 
(|3.45|) has both, charge conjugation C even and odd contributions. The C 
even part corresponds to the "Pomeron" and this is not a Lorentz-vector 
exchange, but the coherent sum of spin 2,4,6,... exchanges (cf. [^7|, |53| ). 
The C-odd part corresponds to the "Odderon" introduced in |S8[] and this is 
indeed a Lorentz vector exchange. 



The quark-quark scattering amplitude ( |3.45 ) is governed by the correla- 
tion function of two connectors or string operators V± associated with two 
light-like Wegner- Wilson lines (Fig. 13). 

The first numerical evaluations of ( |3.45| ) using the methods of the SVM 
were done in |j59"f . However, it turned out that quark-quark scattering was 
calculable from ( |3.45 ) for abelian gluons only. Indeed, we are embarked on a 



program to reproduce in this way the results obtained in |25j in the framework 
of perturbation theory for high energy scattering pD[. For the non-abelian 



gluons difficulties arose having to do with our neglect of quark confinement 
This was really a blessing in disguise and the solution proposed in 



59 



was 



to consider directly hadron-hadron scattering, representing the hadrons as 
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qq and qqq wave packets for mesons and baryons, respectively. We will see 
below how this is done. 

3.5 The scattering of systems of quarks, antiquarks, 
and gluons 

Let us consider now the scattering of systems of partons. As an example we 
study the scattering of two qq pairs on each other: 

q(l) + g(l') + q(2) + q(2') - q(3) + q(3') + g(4) + q(4'), (3.50) 

where we set q(i) = qfa, s u , q(i') = g(p-,s-,AQ (i = 1,...,4) with 
Pi, Si, A4 (p'i, s'i, A'j) the momentum, spin, and colour labels for quarks (an- 
tiquarks). We assume the particles with odd (even) indices i to have very 
large momentum components in positive (negative) x 3 direction (Fig. 14), 
i.e. we assume: 

Pi+,p' i+ — > oo for i odd, 

Pi-,p'i- — > oo for i even. (3.51) 

The transverse momenta are assumed to stay limited. 

Of course, the reduction formula can also be applied for the reaction 
( |3.50|) . We have to be careful in keeping disconnected pieces. The further 



strategy is completely analogous to the one employed in sect. 3.2 for deriving 
(|3.7| )- (|3.17|) . As we dropped the w-channel exchange diagrams in sect. 3.2, 
we drop now all terms which are estimated to give a vanishing contribution 
to high energy small momentum transfer scattering. These terms are char- 
acterized by large momenta of order of the cm. energy yfs flowing through 
gluon lines, leading to suppression factors of order 1/s. Keeping only the 
t-channel exchange terms and performing all the steps as done in sections 
3.2, 3.3 for quark-quark scattering leads finally to a very simple answer for 
the S'-matrix element correponding to reaction ( |3.50|) in the limit ( |3.51| ): 



(3,3',4,4'|S|1,1',2,2'>^ 

([5(3, 1) - iZfM'ufflM*, 10 " iZ^M% lv (G)\ 
[5(4,2) -iZ^M\ 2 (G)][8(A',2') - iZ^M'\, 2 ,(G)]) G . (3.52) 
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Here 5(i,j) is as defined in ( |3.6| ), M^G) and M r A2 (G) are as in ( |3.22|) and 
M%, V (G) and Ai' r 4 , 2 ,(G) are the corresponding amplitudes for the scattering 
of antiquarks on the gluon potential G\(x). We define 



(j'l 



(3.53) 



where v s > } A'.(p'j) is the Dirac and colour spinor for the antiquark q(j'). We 
have then with S r the retarded Green's function for quarks in the gluon 
potential G\(x): 



Ml AG) = -(j'\(i0- m ')S r (i0+m')\k'), 



(3.54) 



where (fc', f) = (3', 1'), (4', 2'). 

In the high energy limit ( 3.51 ) the scattering amplitudes (|3.54 ) can again 
be obtained in the eikonal approximation. Indeed, we can just use C-invari- 
ance to get: 



M'l v (G) -> i^/p' 3+ p' 1+ ■ 5 4 , s ; 
J ds_d 2 xrexp ^(p' 3 — p[)+x- 
[Vl(oo,x-,x. T ) - fl]^,Ai, 



HPs - PiJr • x T 



(3.55) 



(P4 ~ P2)r • yr 

(3.56) 

In section 3.1 we have argued that over the time interval ( |3.3|) we can 
neglect parton production and annihilation processes. For the scattering over 
such a time interval we should have an effective wave function renormalization 
constant 

Zj, = 1, (3.57) 

since the deviation of from 1 is just a measure of the strength of quark 
splitting processes: q —>■ q + G etc. In the calculation of in the framework 
of the SVM to be described below, one finds indeed Z$ — 1, showing the 
consistency of this approach with the simple physical picture of sect. 3.1. 



M'l, 2 ,(G) - 
J dy + d 2 y T exp 
W+(y+, oo,y T ) 



P4-P2- ■ £4,4 
^(P4-^)-yn 
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Anticipating this result we see from Q), ( |3T39| ), ( |3T40D and (|3~55| ), ( g3g ) 
that in the S"-matrix element Q3.52j ) the 5(k,j) (5(k',j')) terms cancel with 
the 11 terms in M r kj (G) (M' r klj ,(G)) in the limit This leads us to the 

following simple rules for obtaining the S"-matrix element in the high energy 
limit: For the right-moving quark (1 — > 3) we have to insert the factor: 

S, + (3,l) = v ^.i,„„/dx^dV 

ex P -p\)+x- - ?(p 3 - Pi)r • x T 

VL(oo,x_,x t )a8^i- ( 3 - 58 ) 
For the right-moving antiquark (1' — > 3') we have to insert the factor: 



Sq + {3', 1') = \/p' 3+ p' 1+ ■ 6^ J dx_d 2 x T 



exp 



-Xp's ~ v'i)+%- ~ *(Ps - Pi)t ■ x T 



V^foo, x_, x t )a' A'- 



(3.59) 



For the left-moving quark (2 — > 4) and antiquark (2' — > 4') we have to 
exchange the + and — labels everywhere in ( |3.58| ) and (|3.59| ). This gives: 



2/T 



exp -(p 4 - V2)-y+ ~ «(P4 - P2)t • Yt 
T4(y + ,oo,y T )A4,A 2 , 



(3.60) 



<V(4', 2') = yy 4 _p' 2 _ ■ 5 si , 4 J dy + d 2 y T 



exp 



;(P4 ~ Ph)-V+ ~ Kp'a ~ p£)t ■ Yt 



^+(y+,oo,y r )x;,^ 



(3.61) 



Finally we have to multiply together the factors S q ± , Sq± and integrate 
over all gluon potentials with the functional integral measure (|3 . 1 4 ) to get 



(3,3',4,4'|S|1,1',2,2'> = (^ + (3,l)5 g - + (3 / ,l / )V(4,2)^--(4',2')) G . (3.62) 



36 



Going from quarks to antiquarks corresponds, of course, just to the change 
from the fundamental representation (3) of SU(3) C to the complex conjugate 
representation (3*) as we see by comparing fl3.58| ) with ( |3.59| ) and ( ^.60[ ) with 



(pi]). 

It is an easy exercise to show that these rules can be generalized in an 
obvious way for the scattering of arbitrary systems of quarks and antiquarks 
on each other. Here we always assume that we have one distinguished col- 
lision axis and that one group of partons moves with momenta approaching 
infinity to the right, the other group to the left. The transverse momenta 
are assumed to stay limited. 

In appendix B we show that these rules can also be extended to gluons 
participating in the scattering. We simply have to change the colour repre- 



sentation in (|3.58|) , ( |3.60|) from the fundamental to the adjoint one. In detail 



we find that for a right-moving gluon 

G(pi,ji,ai) -> G(p 3 , j s ,a 3 ) (3.63) 

the following factor has to be inserted in the S'-matrix element (cf. appendix 
B): 



S G+ (3, 1) = • <Wi 

J dx_d 2 XT exp 

V_(oo,a;_,XT)o3 ) ai- (3-64) 



1 

2 (Ps - Pi)+ X - - *(Ps - Pi)t ■ x T 



Here are the spin indices which are purely transverse, 1 < jj 3 < 2. The 
colour indices are ax,a 3 with 1 < 013 < 8 and V is the connector for the 
adjoint representation of SU(3) C (cf. (|B.23|) ). 



For a left-moving gluon we have again to exchange + and — labels in 

(Pi- 



3.6 The scattering of wave packets of partons 
representing mesons 

In this section we will go from the parton-parton to hadron-hadron scattering. 
Our strategy will be to represent hadrons by wave packets of partons, where 
we make simple "Ansatze" for the wave functions. Then the partonic S'- 
matrix element obtained by the rules derived in section 3.5 will be folded 
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with these wave functions to give the hadronic S'-matrix elements. Of course, 
we always work in the limit of high energies and small momentum transfers. 
Let us start by considering meson-meson scattering: 

M^PO + M 2 (P 2 ) - M 3 (P 3 ) + M 4 (P 4 ), (3.65) 

where are again the right movers, M 2i 4 the left movers. We make simple 
"Ansatze" for the mesons as qq wave packets as follows: 

|M J (P,)) = /dV/dC-^<, s ,(C,P,) 



(2tt) 3 / 2 

P 3 -Pj, 

(j = l,...,4), (3.66) 



-n=S AjtA >.\q{pj, Sj, Aj), q(Pj - Pj , s'j, Aj)) 



where for j = 1,3: 



Pj + — > oo 
Pj+ = 

PjT = oPjT + Pt (3.67) 



and for j = 2, 4: 



Pi- = C^-, 

PjT = 2 P ir + Pt- (3-68) 

Here ( is the longitudinal momentum fraction of the quark in the meson, p<r 
is the relative transverse momentum of q and q. 

We stress that we are not restricting ourselves to spin mesons only. With 
appropriate functions h((, p?) in (|3.66|) we can represent states of gg-mesons 
of arbitrary spin. 

We choose for our states ( |3.66| ) the usual continuum normalization: 

(MjiP'^M^P)) = (27r) 3 2P°5 3 (P' - P). (3.69) 



With ( |3.6|) this requires: 

J d 2 p T J\ 2C(1-C)^(C,Pt)^(C,Pt) = 1 (3.70) 
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In ( |3.69| ) and ( |3.70| ) no summation over j is to be taken. 

For later use we define the wave functions in transverse position space at 
fixed longitudinal momentum fraction (: 

^, s '(C,x T ) := y/2£(l - C)-^ j d 2 p T exp(ip T • x T )^ s ,(C, p T ) 

(7 = 1,. ..,4). (3.71) 

With this we define profile functions for the transitions Mj — > for right 
and left movers as: 

^(x T ) := / 1 dC« s ,(C,x r ))V^(C,x T ), (3.72) 

J 

where k,j are both odd or even. Clearly we have (cf. ( |3.70|) ) 

J d 2 x T w j:j (x T ) = 1, 

(no summation over j). (3.73) 

Let us first study the transition of the right movers alone, i.e. the "reac- 
tion" : 

M 1 {P l ) -> M 3 (P 3 ). (3.74) 

For stable mesons Mi i3 we should find that the corresponding S'-matrix ele- 
ments are identical to the matrix elements of the unit operator. Is this borne 
out in our approach? 

From the rules given in sect. 3.5 we find easily the S'-matrix element for 
the transition 

g(l) + g(l')-g(3)+g(3') (3.75) 

in the form 

(3,3'|S|1,1'> = (S 9+ (3,1)%(3',1')> . (3.76) 
After folding ( |3.76| ) with the mesonic wave functions ( |3.66| ) we get: 
(M 3 (P 3 )|S|M 1 (P 1 )) = (2vr) 3 2P 1 °5 3 (P 3 - P x ) 

[d 2 z T [ d( 3 [ dCi ^,V(C3,Zt)^ jS /(Ci,Zt) 
J Jo Jo 

(iTr[V_(oo, z-, z T )Vl(oo, 0, 0)]) G . (3.77) 
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Now we remember that we consider the limit P\ + — > oo. Therefore we 
perform a change of variables in the Z- integral by setting: 



This gives: 



z'_ := \Pi + Z-. (3.78) 



(M 3 (P 3 )|5|M 1 (P 1 )) = (2vr) 3 2P 1 ° ( 5 3 (P3 - Pi) 

/ d 2 z T [ d( 3 [ dCi ^.VCCs.Zr)^ ,(Cl,Zr) 
J jo Jo 

/^exp[z(C 3 -CiR] 
(\Tr[V-(oo,-?-z'_,z T )Vl(oo,0,0)}) G 

— > (27r) 3 2P 1 °5 3 (P 3 - Pi) | dV u; 3 ,i(zt) 
(iTr[F„(oo,0,z T )K!(oo,0,0)]) G for P 1+ -> oo. (3.79) 



In ( |3.79|) V_(oo, 0, zr)(V r i(oo, 0, 0)) is the quark (antiquark) connector 



taken along the line C q (Cq), where: 

C q : r -> z 9 (t) = | z T | , (-oo < r < oo), (3.80) 




Cg-.r^ Zq(r) = , (-oo < r < oo). (3.81) 



But from sect. 2.1 we know that the connector Vl along Cq is equal to 
the connector V- taken along the oppositely oriented line Cq (cf. Q2.25|) ). 
Now we will allow ourselves to join the lines C q and Cq at large positive 
and negative times. We can imagine the gluon potentials to be turned off 
adiabatically there. We obtain then from the product of the connectors in 
(|3.79|) a connector taken along a closed lightlike Wegner- Wilson loop 

^Tr[V_(oo, 0, z T )V.!(oo, 0, 0)] — > W+(^z T , z T ). (3.82) 
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Here we define 

W ± (y r ,z r ) = \TrP exp[-ig [ dx^(x)] (3.83) 

o J C± 

with C + (CJ) a lightlike Wegner- Wilson loop in the plane x_ = 0(x + = 0), 
where in the transverse space the centre of the loop is at yx and the vector 
from the antiquark to the quark line is zt (Fig. 15). 

Inserting now everything in ( 3.79|) we get the simple answer: 

(M 3 {P 3 )\S\M 1 (P 1 )) = (2 7 r) 3 2P 1 5 3 (P 3 - Px) 

J d 2 z T w 3A (z T )(W + (^z T ,z T )) G . (3.84) 

In the next section we will evaluate the functional integral in (|3.84j ) in the 
SVM. We will find 

(W + (y T ,z T )) G \ SVM = l. (3.85) 
Inserting this in ( |3.84j ) leads to the expected result (cf. (|3.66j) - (|3.72| )): 

(M 3 (P 3 )|5|M 1 (P 1 )) = (M 3 (P 3 )|M 1 (P 1 )). (3.86) 

In our approach the qq pair in the right-moving meson Mi — > M 3 does 
not interact. Of course this is only valid over our finite time interval (3.3)! 

The techniques developed thus far are now easily employed for the reac- 
tion (|3.65|) . After performing similar steps as above we arrive at the following 
S'-matrix element: 

S fi = 5 fl + t(2'K) i 5(P 3 + P 4 -P 1 -P 2 )T fl , 
T fi = (M 3 (P 3 ),M 4 (P 4 )|T|M 1 (P 1 ),M 2 (P 2 )) 

= -2is J d 2 6 T d 2 £ T dVe iqT ' bT w 3i i(x T )w 4 ,2(yT) 

(W + (ib T ,x T )W_(-^b T ,y T ) - 1) G . (3.87) 

Here s = (Pi + P 2 ) 2 is the cm. energy squared and qr is the momentum 
transfer, which is purely transverse in the high energy limit: 

q T = (Pi - P 3 ) T . (3.88) 
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From ( |3.87|) we see that the amplitude for soft meson-meson scattering at 
high energies is governed by the correlation function of two lightlike Wegner- 
Wilson loops, where one is in the hyperplane x_ = 0, the other in x + = 
0. The transverse separation between the centres of the two loops is given 
by by, the impact parameter. The vectors and give the extensions 
and orientations of the loops in transverse space (Fig. 16). The loop-loop 
correlation function has to be integrated over all orientations of the loops in 
transverse space with the (transition) profile functions of the mesons, u^i 
and u>4,2- Finally a Fourier transform in the impact parameter has to be 
done. 

The methods presented here for meson-meson scattering can of course 
also be employed for scattering reactions involving baryons and antibaryons. 
This is sketched in appendix C. 

3.7 The evaluation of scattering amplitudes in the Minkowskian 
version of the stochastic vacuum model 

In the previous section we have derived expressions for the amplitudes of 
soft meson-meson scattering at high energies in terms of correlation func- 
tions of light-like Wegner- Wilson loops. The task is now to evaluate the 
corresponding functional integral ()a in ( |3.87|) . Surely we do not want to 



make a perturbative expansion there, remembering our argument of sect. 1 



(cf. ( 1.14 )). Instead, we will turn to the SVM which did quite well in its 



applications in Euclidean QCD (sect. 2). Of course, the generalization of the 



SVM to Minkowski space-time is a bold step which was done in ||59 |, [45 
The authors of these refs. proposed to use in Minkowski space-time just the 
assumptions 1-3 of the SVM (cf. sect. 2.4), but after having made a suit- 
able analytic continuation. In this way we obtain for instance the Minkowski 
version of Ass. 1 as (cf. (|2^q) , fl2~HD): 



Ass. 1: The correlator of two field strengths, shifted to a common 
reference point y, is independent of the connecting path and given by: 

(£^G; u (y,x-,C x )G b pa (y,x'-,C x/ )) G = ^F^z), 

Fpupaiz) = ^G 2 {(g^ p g ua - g tl(T 9up)^D(z 2 ) 
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Id d 
+-\— (z a g w - z p g pa ) + — (z p g ua - z a g vp ) ] ■ (1 - k)Di(z 2 )}. 



2 L dz" 



(3.89) 



Here z = x — x' and D(z 2 ) and Di(z 2 ) are denned as in (|2.65 ), ( |2.67 ) for 
z 2 < and by analytic continuation for z 2 > 0. 

As a first application let us calculate the expectation value of one lightlike 
Wegner- Wilson loop. We have from ( 3.83 ) using the non-abelian Stokes 
theorem: 

(W + (y T ,z T )) G = (^TVPexp [-if? / dx p G^x)}) G 

9(x M , a; 



(— Tr Pexp [ — i— 



2 Js+ 



dudv - 



d(u, v) 



G^(R,x;C x )]) G . (3.90) 



Here S + is the (planar) surface spanned into C + (Fig. 15) and parametrized 
by 



x»(u, v ) = u n\ + - (v - i)z M , 

— oo < w < oo, < v < 1, 
where n + is as in ( 3.42j ) and 



(3.91) 



(3.92) 



The reference point on the surface S+ is denoted by R and C x are straight 
lines running from x to R. From ( |3.91|) we find 





( 


) 


















\ 


o J 







d{x> x : x v ) 
d(u, v) 



z^n" - ntz v . 



(3.93) 



Now we apply the cumulant expansion formulae (cf. sect. 2.3) to ( |3.9C| ) 
and use assumptions 1-3 of the SVM. This is completely analogous to the 
calculations done in sect. 2.5. We get: 



(W + (y T , zt))g = exp{ — ^— dudv du'dv'K 2 {x - x' 

6 Js+ Js+ 



+ higher cumulant terms}, 



(3.94) 
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where 



Ko(x 



X 



d(x»,x u )d(x' p ,x' a ) 
d(u,v) d(u',v') 



X 



X{U, V 



X 



III l\ 

x (u , V } 



(3.95) 



It is an easy exercise to evaluate K 2 (x — x') using ( |3.93| ) and F MJ , pCT from 



(137891) . The result is 



K 2 (x-x')=0 (3.96) 
for all x,x' G S+. With assumption 3 of the SVM (cf. (^J^)-(ME)) a11 



higher cumulants are related to the second one. Thus, ( |3.96| ) implies also the 
vanishing of all higher cumulant terms in ( |3.94j ) and we get 



(W+(y T ,z T )) G -- 
A similar argument leads, of course, also to 

(>V_(yT,z T )) G = 



(3.97) 



(3.98) 



In sect. 3.6 we have already used the result ( |3.97| ) in the discussion of the 
transition Mi — > M 3 to obtain ( |3.86|) . We can also use it for the calculation of 
the wave function renormalization constant Z^. The expression we obtained 
for Zrf in (|3.48|) can be interpreted as the expectation value of the non- 
abelian phase factor picked up by a very fast right-moving quark. Now, 
isolated quarks do not exist. The best approximation for it we can think of 
is a fast, right-moving quark-antiquark pair with the antiquark being very 
far away from the quark in transverse direction. In this way we obtain for 
Z^p instead of (|3.48|) : 



lim (-Tr[F_(oo,0,0)yi(oo,0,z T )]) 

Zy|— >oo 3 

lim (W + (--z T , -z T )) G 

1. 



G 



(3.99) 



This result was already used in sect. 3.5, ( p.57p ff. 
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We come now to the evaluation of the loop-loop correlation function of 
(PZD: 



(W + (ib r ,x T )W_(-ib r , y r ) - 1) G 

= ([W+(^b T ,x T ) - l][W_(-ibr,y r ) - 1])g- (3.100) 



Here we used ( p.97| ), ( |3.98| ). The strategy is as before. We want to transform 



the line integrals of W± (cf. (|3.83|) ) into surface integrals using the non 
abelian Stokes theorem of sect. 2.2. Following the authors of [|59j] , [i5 we 
choose as surface with boundary C + and C_ a double pyramid with apex 
at the mid-point of C + and C_ which is the origin of our coordinate system 
(Fig. 17). The mantle of this pyramid is V + + V- and we have with suitable 
orientation 

d(V + + V-) = C+ + C_. (3.101) 

In the transverse projection of Fig. 16 the basis surface S + (S-) appears 
as the line q + q + (?-?_) and the mantle surface V + (VJ) as the triangle 
0g + g + (0g_g_). From the non-abelian Stokes theorem we obtain now for 
(1T00D: 

([W + (~b r ,x T ) - l][W_(-ibr,y T ) - 1]) G 



({lTrPexp[-4 / dMd /(^,x-) ^ _ 
3 2 Jr+ oiu.v) 



d(u, v) 

{^TrPexp[-zf / du'dv' d{ f; ^ G pa (0, x'; C*)] - 1}) G . 
3 2 o{u,v') 

(3.102) 

So far, it has not been possible to use some version of the cumulant expansion 
for (|3.102|) . Thus in [[59], |45| the path-ordered exponentials on the r.h.s. of 



( |3.102[ ) were expanded directly. The structure of this expansion is as follows: 



([W + -1][W_-1]> 



G 



3 JT+ 



V+ JV+ 



GG 



(3.103) 
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The trace of a single shifted field strength vanishes. This means that we 
cannot exchange a single coloured object in meson- meson scattering ( 3.65|) . 



The first non-trivial contribution in ( |3.103| ) comes from the term with two 
shifted field strengths in each trace. The corresponding correlation function 

I GG\ ~Tr[/ / GG\) G 
3 Jv+ Jv+ 3 Jv- JV- 

can now be evaluated using the assumptions 1-3 of the Minkowskian ver- 
sion of the SVM. In doing so it is advantageous to transform the surface 
integrals over the mantles of the pyramids V± into surface integrals over S± 
and integrals over the volumes V± enclosed by V+ and S+ and V- and S-, 
respectively: 

d(V±) =P±-S±. (3.104) 

Then one can use the ordinary Gauss theorem to get simpler integrals. 

We will not enter into the details of these calculations here, but only note 
that the integrations along the directions x + and x_ can easily be done ana- 
lytically and that one ends up with integrals over the projections of S± and V± 
into the transverse space. Thus, one finally needs the correlator (|3.89f) for 
space-like separations only: 



[x 



X 



'f = z 2 


<o, 






) 




( 





XT , 


x' = 






o / 




\ 






(3.105) 



where x T runs over the triangle A + = 0q + q + and x' T over A_ = 0q_q_ in 
Fig. 16. For space-like separations the correlator functions D(z 2 ) and Di(z 2 ) 
in Q3.89p are as in Euclidean space time. The resulting expressions are then 



of the following structure 

([W + -1][W_-1]) G ~ 

{ / d 2 z T [ d 2 z' T G 2 [..D(-(z T - z' T ) 2 ) + ...^(-(zt - z^) 2 )]} 2 . 

(3.106) 

These integrals have to be evaluated numerically. 
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With the methods outlined above, we have obtained an (approximate) 
expression ( |3.106 ) for the functional integral ()g governing the meson-meson 
scattering amplitude ( |3.87| ). Note that the nonperturbative gluon condensate 
parameter Gi sets the scale in ( |3.106|) and in the integrals to be performed 
there the vacuum correlation length a enters through the D and D\ functions. 
Thus, on dimensional grounds, we must have in our approximation: 



<(W+ - 1)(W_ - 1)> 



= <W!* 5:, *£ 

\ a a a , 



(3.107) 



where / is a dimensionless function. To obtain the meson-meson scattering 
amplitude Q3.87D we still have to integrate over the profile functions w 3j i(x T ) 
and ^4,2 (yr)- Here the transverse extensions of the mesons - i.e. of the wave 
packets representing them - enter in the results. 

For a detailed exposition of the numerical results obtained in the way 
sketched above we refer to [|45|| . Here we only discuss the outcome for proton- 
proton scattering when treated in a similar way (cf. appendix C). A fit to the 
numerical results gives for the total cross section and the slope parameter at 
t = of elastic proton-proton scattering the following representation: 



<7tot(PP) 

d 



0.00881 [ 

a 



3.277 



(3tt 2 G 2 



,10 



(3.108) 



J vv 



1.558a 2 + 0.454i?;. 



Here R p is the proton radius and the formulae ( |3.108| 

1 < Rp/a < 3. 



(3.109) 
( |3.109| ) are valid for 
(3.110) 



To compare (|3.108|) , ( p.!09| ) with experimental results, we can, for in- 
stance, consider the cm. energy y/s = 20 GeV and take as input the follow- 
ing measured values (cf. p5|): 



Pomeron part 
b 



PP ~ 

R p = 



35 mb, 
12.5 GeV" 2 , 
Rp,elm = 0.86 fm. 



We obtain then from flSTTDgp and (|3~T09|) : 

a = 0.31 fm, 
G 2 = (507 MeV) 4 . 



(3.111) 



(3.112) 
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The values for the correlation length a and for the gluon condensate G2 come 
out in surprising good agreement with the determination of these quantities 
from the fit to the lattice results ( |2.70| ). 

But perhaps we were lucky in picking out the right cm. energy ^/s 
and radius for our comparison of theory and experiment. What about the 
s-dependence of the total cross section a to t and slope parameter 6? The 
vacuum parameters G 2 and a should be independent of the energy yfs. On 
the other hand, from the discussion in sect. 3.1 leading to (|3~2|), it seems 
quite plausible to us that the effective strong interaction radii R of hadrons 
may depend on y^. Let us consider again pp (or pp) elastic scattering. Once 
we have fixed G2 and a from the data at y/s = 20GeV (|3.108|) and ( |3.109| ) 
give us (Jtotipp) and b pp in terms of the single parameter R p , i.e. we obtain 
as prediction of the model a curve in the plane b pp versus cr tot (pp). This is 
shown in Fig. 18. It is quite remarkable that the data from -^i = 20GeV up 
to Tevatron energies, ^/s = 1.8 TeV follow this curve. 

Summarizing this section we can say that explicit calculations for high 
energy-elastic hadron-hadron scattering near the forward direction have been 
performed combining the field-theoretic methods of |53| and the Minkowski 
version of the stochastic vacuum model of [59fl,[|4l)|. The results are encourag- 
ing and support the idea that the vacuum structure of QCD plays an essential 
role in soft high-energy scattering. We want to point out that these calcu- 
lations also resolve a possible paradoxon of QCD: On the one hand there 
are suggestions that the gluon propagator must be highly singular (probably 
oc (Q 2 )~ 2 , cf. e.g. [31]) for momentum transfers Q 2 — > in order to pro- 
duce confinement. On the other hand high energy scattering amplitudes are 
completely regular for t — > 0. A singular gluon propagator will lead in the 
2-gluon exchange model to a singularity for t = not only for quark-quark 
scattering but also for hadron-hadron scattering if the latter are considered 
e.g. as colour dipoles. The resolution of the paradoxon which we can present 
is intimately connected with the confinement mechanism which we found in 
the SVM in Sect. 2. The short range correlation of the gluon field strengths 
governs the t-dependence of the hadronic scattering amplitudes and gives rise 
to their regularity for t = 0. The gluon propagator on the other hand can be 
singular for Q 2 — > 0, since the gluon potentials have a long range correlation 
as we have seen in Sect. 2. The result (|3.108| ) for the total cross section de- 
pends also on the proton radius. This radius dependence does not saturate 
for large radii in the calculation with non-abelian gluons but does saturate 
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in an abelian model [|l|. Thus with non-abelian gluons we do not get the 
additive quark model result [17], and thus not the picture of Donnachie and 



Landshoff |Z3|| , where the "soft Pomeron" couples to individual quarks in the 
hadrons. The strong radius dependence in ( |3.108|) is due to the D-term in the 
correlator (|3.89Q which is related to the effective chromomagnetic monopole 
condensate in the QCD vacuum and which gives rise to the linearly rising 
quark- ant iquark potential, i.e. to string formation (see sect. 2.5). The calcu- 
lations reported above suggest that in high energy scattering this same term 
gives rise to a string-string interaction which leads to the radius dependence 
in p. 108|) . Note that one does not have to put in the strings by hand. They 
enter the picture automatically through our lightlike Wegner-Wilson loops. 
The radius dependence of the cross section occurs, of course, also for meson- 
baryon and meson-meson scattering and gives a quantitative understanding 
of the difference between the Kp and 7rp total cross sections and slope param- 
eters at high energies. For this and for further results we refer to [fj5] |, [p2^] . 
The success of the calculation for pp{pp) scattering describing correctly the 
relation of the total cross section versus the slope parameter from a/s ~ 20 
GeV up to y/s = 1.8 TeV suggests the following simple interpretation: In soft 
elastic scattering the hadrons act like effective "colour dipoles" with a radius 
increasing with cm. energy. The dipole-dipole interaction is governed by 
the correlation function of two lightlike Wegner-Wilson loops which receives 
the dominant contribution from the same non-perturbative phenomenon - 
effective chromomagnetic monopole condensation - which leads to string for- 
mation and confinement. 
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4 "Synchrotron Radiation" from the 

Vacuum, Electromagnetic Form Factors 
of Hadrons, and Spin Correlations in the 
Drell-Yan Reaction 

Let us consider for deflniteness again a proton-proton collision at high cm. 
energy i/i ^> m p . We look at this collision in the cm. system and choose as 
£ 3 -axis the collision axis (Fig. 19). According to Feynman's parton dogma 
|63|j the hadrons look like jets of almost non-interacting partons, i.e. quarks 
and gluons. Accepting our previous views of the QCD vacuum (Sect. 2), 
these partons travel in a background chromomagnetic field. 

What sort of new effects might we expect to occur in this situation? 
Consider for instance a quark- ant iquark collision in a chromomagnetic field. 
In our picture this is very similar to an electron-positron collision in a storage 
ring (Fig. 20). We know that in a storage ring e~ and e + are deflected 
and emit synchrotron radiation. They also get a transverse polarization due 
to emission of spin-flip synchrotron radiation [Ell, |65||. Quite similarly we 



can expect the quark and antiquark to be deflected by the vacuum fields. 
Since quarks have electric and colour charge, they should then emit both 
photon and gluon "synchrotron radiation" . Of course, as long as we have 
quarks within a single, isolated proton (or other hadron) travelling through 
the vacuum no emission of photons can occur, and we should consider such 
processes as contributing to the cloud of quasi-real photons surrounding a 
fast-moving proton. (This is similar in spirit to the well-known Weizsacker- 
Williams approximation.) But in a collision process the parent quark or 
antiquark will be scattered away and the photons of the cloud can become 
real, manifesting themselves as prompt photons in hadron-hadron collisions . 



In ref. we have given an estimate for the rate and the spectrum of 
such prompt photons using the classical formulae for synchrotron radiation 
[p5[| . A more detailed study of soft photon production in hadronic collisions 



was made in pB |. A sketch of our arguments and calculations is as follows. 

In Sect. 2 we discussed the domain picture of the QCD vacuum. In 
Euclidean space time we have domains in the vacuum of linear size ~ a. 
Inside one domain the colour fields are highly correlated. The colour field 
orientations and domain sizes fluctuate, i.e. have a certain distribution. If 
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we translate this picture naively to Minkowski space, we arrive at colour 
correlations there being characterized by invariant distances of order a. Then 
the colour fields at the origin of Minkowski space, for instance, should be 
highly correlated with the fields in the region 

|x 2 |<a 2 (4.1) 

(cf. Fig. 21). Consider now a fast hadron passing by with one of its quarks 
going right through x = on a nearly lightlike world line. It is clear that 
in such a situation the quark will, from the point of view of the observer, 
spend a long time in a correlated colour background field. An easy exercise 
shows furthermore that two quarks of the same hadron will have a very 
good chance to travel in two different colour domains. The argument is in 
essence as follows. The quarks have a transverse separation of the order of 
the hadron radius R whereas the transverse size of a domain is of order a 



and we have a /R <C 1 (cf. (|2.72|) ). Each quark will then wiggle around 
due to the deflection by the background colour fields in an uncorrelated 
fashion. This gives us a justification for adding the synchrotron gluon and 
photon emission of the quarks incoherently . The result we found can be 
summarized as follows: In the overall cm. system of the hadron-hadron 
collision "synchrotron" photons should appear with energies uo < 300 — 500 
MeV, i.e. in the very central region of the rapidity space. The number 
of photons per collision and their spectrum are — apart from logarithms 
- independent of the cm. energy y/s. The dependence of the number of 
photons on the energy uj and on the emission angle •&* with respect to the 
beam axis is obtained as follows for pN collisions: 

d ^ = SM 2/3 -S(cosr) (4.2) 



dcudcos'i?* u 1 / 3 

Here a is the fine structure constant and Z e ff is the length or time over which 
the fast quark travelling in the background chromomagnetic field B c obtains 
by its deflection a transverse momentum of order px ~ 300 MeV, the mean 
transverse momentum of quarks in a hadron (Fig. 22): 

Z« = ^r- (4.3) 
gB c 



The quantity S in (|4.2| ) sums up the contributions from all quarks of the 



initial and final state hadrons. It involves an integration over the quark 
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distribution functions of these hadrons. In [RH1 we found 



for pAT collisions at y/s = 29 GeV. 



Our result (|4.2|) for synchrotron photons should be compared to the inner- 

« — ^ (4.5) 



bremsstrahlungs spectrum 

dn. 



7 



a; sin $* 

bremsstr. 



dcvd cos $* 

The "synchrotron" radiation from the quarks ( |4.2| ) is thus harder than the 
hadronic bremsstrahlung spectrum. This is welcome, since for uo — > brems- 
strahlung should dominate according to Low's theorem JB7| . 

It is amusing to note that in several experiments an excess of soft prompt 
photons over the bremsstrahlung calculation has been observed [)58[]-||72jl. The 
gross features and the order of magnitude of this signal make it a candidate 
for our "synchrotron" process. A detailed comparison with our formulae has 
been made in pB| for the results from one experiment [72] with encouraging 



results. This is shown in Fig. 23 for the kr spectrum of photons at y = 0, 
where 

= u>sini?*, 

y = -lntaa(i?*/2). (4.6) 

We see that the addition of synchrotron photons to the bremsstrahlung ones 
improves the agreement of theory with the data considerably. We deduce 
from Fig. 23 l cS ~ 20 - 40 fm. Taking l eS = 20 fm and p T = 300 MeV, we 
find for the effective chromomagnetic deflection field from QL3|) 



gB c =^ = (55 MeV) 2 . (4.7) 

left 



This is much smaller than the vacuum field strength (2.7). Our interpretation 
of this puzzle is as follows: The colour fields in a fast moving hadron must be 
shielded. Indeed, a chromomagnetic field of the strength ( |2.7| ) would lead to 
a ridiculously small value for the radius of cyclotron motion of a fast quark. 
The necessary shielding could be done by gluons in a fast hadron. We know 
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from the deep inelastic lepton-nucleon scattering results that a fast nucleon 
contains many gluons. We may even be brave and turn the argument around: 
in order for a fast hadron to be able to move through the QCD vacuum, the 
very strong vacuum chromomagnetic fields must be shielded, making gluons 
in a fast hadron a necessity . Thus soft photon production in pp collisions 
may give us a quite unexpected insight into the quark and gluon structure 
of fast hadrons. 

The next topic we want to discuss briefly concerns electromagnetic form 
factors of hadrons. 

We have argued above that the colour fields in the vacuum should give 
a contribution to the virtual photon cloud of hadrons and we made an es- 
timate of the distribution of these photons using the synchrotron radiation 
formulae. Consider now any reaction where a quasi- real photon is emitted 
from a hadron with the hadron staying intact and the photon interacting 
subsequently. In Fig. 24 we draw the corresponding diagram for a nucleon 
N: 

N(p)-,N(p') +1 (q). (4.8) 
The flux of these quasi-real photons is well known. The first calculations 



in this context are due to Fermi, Weizsacker, and Williams JF3). For us the 
relevant formula is given in eq. (D.4) of J73J]. Let E be the energy of the initial 
nucleon, G^(Q 2 ) its electric Sachs formfactor, and let u and q 2 = —Q 2 be the 
energy and mass of the virtual photon. Then the distribution of quasi-real 
photons in the fast-moving nucleon is given by 



dn^ l) = --^\G N E (Q 2 )Y (4.9) 
where we neglect terms of order oj/E and Q 2 /m 2 N and assume 

Q 2 » QLn - (4.io) 

We call (|4.9| ) the exclusive flux since the nucleon stays intact. Now we want 
to translate ( [4.9| ) into a distribution in u and the angle of emission of the 
7 (cf. Fig. 24). A simple calculation gives 

Q 2 ~ u 2 sm 2 ^\ (4.11) 



2 



2 
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dni excl) 



2a du d$* 



G%(u 2 sm 2 $*) 



(4.12) 



7r to sint?* 

Now we made an "exclusive-inclusive connection" argument in [j(J: We re- 
quire dn( excl ) /da;d cos i?* to behave as u;~ 1//3 for fixed $* as we found in (|4.2|) . 
This implies for the form factor G E (Q 2 ) a behaviour as (Q 2 ) 1 ^ 6 - 

Thus we arrive at the following conclusion: The proton form factor 
G P E should contain in addition to a "normal" piece connected with the to- 
tal charge and the hadronic bremsstrahlung in inelastic collisions a piece 
oc (Q 2 ) 1 / 6 connected with "synchrotron" radiation from the QCD vacuum. 
For the neutron which has total charge zero we would expect the "normal" 
piece in G E to be quite small and the "anomalous" piece to be quite impor- 
tant for not too large Q 2 . Thus the neutron electric formfactor should be an 
interesting quantity to look for "anomalous" effects oc (Q 2 ) 1 ^ 6 . 

In Fig. 25 we show the data on the electric formfactor of the neutron from 



We superimpose the curve 



G" syn) (Q 2 ) = 3.6-10- 



1/6 



5 fm 



-2 



(4.13) 



which is normalized to the data at Q 2 = 5fm -2 . We see that except in the 
very low Q 2 region we get a decent description of the data. For Q 2 — > ( |4.13j ) 
has to break down since G E (Q 2 ) is regular at Q 2 = 0. Indeed one knows the 
slope of G E (Q 2 ) for Q 2 = from the scattering of thermal neutrons on 
electrons ( |f77f and references cited therein): 



dGUQ 2 ) 



dQ 2 



0.019 fm" 



(4.14) 



We see from Fig. 25 that the behaviour of G^(Q 2 ) has to change rather 
quickly as we go away from Q 2 = 0. We will now make a simple ansatz 
which takes into account that G^(Q 2 ) can have singularities in the complex 
Q 2 -plane only for 

-oo < Q 2 < -4m 2 , 

where m n is the pion mass. We require a (Q 2 ) 1 ^ 6 behaviour for positive Q 2 
and take the slope of C%{Q 2 ) at Q 2 = from experiment Q4.14| ). This leads 
us to the following functional form for G^(Q 2 ): 

-5/6 



G n E (Q 2 ) = 0.019 fm 2 ■ Q 2 



1 + 



4 m 2 



(4.15) 
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It is amusing to see that this gives a decent description of the data (Fig. 25). 

What about the electric form factor of the proton G P E (Q 2 )1 Here, clearly, 
we have a dominant "normal" piece connected with the total charge. We 
will assume that this normal contribution is represented by the usual dipole 
formula 



G D (Q 2 



1 



Q 2 



m 



d , 



m 



D 



18.23 fm~ 2 =0.710 GeV 2 



(4.16) 



which gives a good representation of the data for Q 2 = 2—4 GeV 2 J78"| . Let us 
add to this an anomalous piece for smaller Q 2 , connected with synchrotron 
radiation, and let us assume that this is a purely isovector contribution, 
consistent with the singularity at Q 2 = —4m 2 in ( |4.15|) . We obtain then 
from ( [4.15| ) the following ansatz for G E (Q 2 ): 



G P E (Q 2 ) = G D (Q 2 ) - G E (Q 2 ) = G D (Q 2 )(1 - A(Q 2 )), 



where 



A(Q 2 ) = 0.019 fm 2 ■ Q 2 (l + 



Am 2 



-5/6 



Q 2 



rn 



(4.17) 



(4.18) 



D 



We predict a deviation of the ratio G P E (Q 2 ) j G d{Q 2 ) from unity for small 
Q 2 . It is again amusing to note that such a deviation is indeed observed 
experimentally j|7, 48| . Our ansatz does even quite well quantitatively (Fig. 
26). For the electromagnetic radius of the proton we predict from ( [4.17] ) 



-6 
12 



dG p E (Q 2 ) , 
dQ 2 Iq 2 =o 

+ 6 ■ 0.019 fm 2 



D 



= (0.88 fm) 2 . 

This checks well with the experimental values quoted in H7 



(4.19) 



(4) 1/2 = 

depending on the fit used for G P E (Q 2 ) at low Q 2 



0.88 ± 0.03 fm, or 
0.92 ± 0.03 fm, 



(4.20) 
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To conclude this brief discussion of nucleon form factors, we can summa- 
rize our picture as follows: The quarks in the nucleon make a cyclotron-type 
motion in the chromomagnetic vacuum field. This leads to a spreading out 
of the charge distribution of the neutron 

p„(x) oc |x|- 10 / 3 (4.21) 

corresponding to G^(Q 2 ) oc (Q 2 ) 1 ^ 6 . The same effect should lead to a devia- 
tion of the proton form factor G P E (Q 2 ) from the dipole formula for small Q 2 . 
Concerning the sign of G^Q 2 ) we can in essence follow the arguments put 



forward in |79 



One might think — maybe rightly — that these ideas are a little crazy. 
But we have also worked out some consequences of them for the Drell-Yan 
reaction ( |1.9|) , which make us optimistic. In the lowest order parton process 
contributing there, we have a quark-antiquark annihilation giving a virtual 
photon 7*, which decays then into a lepton pair (Fig. 1): 

q + q -> 7* -> tt. (4.22) 

In the usual theoretical framework q and q are assumed to be uncorrelated 
and unpolarized in spin and colour if the original hadrons are unpolarized. 
From our point of view we expect a different situation. Let us go back to Fig. 
21 where we sketched the world line of a quark of one fast hadron in a colour 
domain of extension \x 2 \ ^ a 2 . Let the quark q and antiquark q in ( |4.22j ) 
annihilate at the point x = 0. Here q and q come from two different hadrons 
hi, h 2 moving with nearly light-like velocity in opposite directions. It is clear 
that in this situation q and q will spend a long time in a highly correlated 
colour background field (Fig. 27). In |l4j] we speculated that this may lead to 
a correlated transverse spin and colour spin polarization of q and q due to the 



chromomagnetic Sokolov-Ternov effect |64], |65|. In []80[ we worked out this 
idea in more detail and found that a transverse qq spin correlation influences 
the £ + £~ angular distribution in a profound way. Then our colleague H. J. 
Pirner pointed out to us that data which may be relevant in this connection 
existed already And very obligingly these data show a large deviation 



from the standard perturbative QCD prediction. On the other hand, we can 
nicely understand the data in terms of our spin correlations and thus vacuum 
effects in high energy collisions. For more details we refer to [|80f| . If such 
spin correlations are confirmed by experiments at higher energies, we would 
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presumably have to reconsider the fundamental factorization hypothesis for 
hard reactions which we sketched in Sect. 1 and which is discussed in detail 
in PL 



5 Conclusions 

In these lectures we have discussed various ideas connected with non-pertur- 
bative QCD and in particular with the QCD vacuum structure. In Sect. 2 
we introduced connectors, the non-abelian Stokes theorem and the cumulant 
expansion. Then we presented the domain picture of the QCD vacuum and 
the stochastic vacuum model (SVM). The latter is consistent with the view 
of the QCD vacuum acting like a dual superconductor: We found that in the 
SVM the vacuum contains an effective chromomagnetic monopole conden- 
sate, whose effect is parametrized by the D-term in the gluon field strength 
correlator (|2.61|) . With these tools we could calculate the expectation value 
of the Wegner- Wilson loop in the SVM. We found a linearly rising potential 
between a heavy quark- ant iquark pair, QQ. This is related to the formation 
of a "string" , a chromoelectric flux tube between QQ as can be seen explic- 
itly in the SVM |52]| . Thus in this framework confinement is an effect of the 
nontrivial QCD vacuum structure. All calculations in sect. 2 were done in 
Euclidean space time. 

In appendix A we discuss some problems which arise when one considers 
higher cumulant terms in the SVM. We propose as remedy for these problems 
to do the calculation of the expectation value of the Wegner- Wilson loop in 
an iterative way, summing step by step over the contributions of various 
plaquettes. This also leads us to a proposal for including dynamical fermions 
in the calculation with the expected result that the linearly rising potential 
levels off and goes to a constant at large QQ separation, where, of course, 
we have then mesons Qq, qQ formed of a heavy quark Q (antiquark Q) and 
a light antiquark q (quark q). 

In sections 3, 4 we looked at various possible effects of the nontrivial 
vacuum structure in the Minkowski world. In section 3 we gave a detailed 
account of a field-theoretic method for the calculation of scattering ampli- 
tudes of high-energy soft hadronic reactions. We started from the functional 
integral and made high energy approximations in the integrand. This led us 
to give general rules for writing down scattering amplitudes for high energy 
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soft reactions at the parton level. Then we considered as representation for 
hadrons wave packets of partons. The corresponding scattering amplitudes 
for (quasi-)elastic hadron-hadron scattering were found to be governed by 
the correlation functions of lightlike Wegner- Wilson loops. The evaluation of 
these correlation functions was possible with the help of the Minkowski ver- 
sion of the stochastic vacuum model. The comparison with experiment gave 
very good consistency, supporting the view that the QCD vacuum struc- 
ture plays an essential role in high energy soft reactions. The framework 
developed in these lectures can be applied directly to elastic and diffractive 
hadron-hadron scattering at high energies. In principle we should also be 
able to apply it to non- diffractive reactions, fragmentation processes etc., 
but this remains to be worked out. 

In sect. 4 we argued that some more startling QCD vacuum effects in 
high energy collisions may be the appearance of anomalous soft photons in 
hadron-hadron collisions due to "synchrotron radiation" and spin correlations 
in the Drell-Yan reaction due to the chromodynamic Sokolov-Ternov effect. 
Furthermore, we gave an argument that electromagnetic form factors at small 
Q 2 should reflect the vacuum structure. Finally we would like to mention 



that in [82] the rapidity gap events observed at HERA are quantitatively 
described in terms of the parton model but invoking again nonperturbative 
QCD effects, possibly connected with the vacuum structure. Another place 
where the QCD vacuum structure may show up is in certain correlations 
of hadrons in Z° decays to two jets fl4|, [83|] for which there is also some 



experimental evidence [Q . 

We hope to have convinced the reader that the non-perturbative structure 
of the QCD vacuum is useful in order to understand confinement of heavy 
quarks. In our view this vacuum structure manifests itself also in high energy 
soft and hard reactions. We think it is very worth-while to study such effects 
both from the theoretical and the experimental point of view. 
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Appendix A: 

Higher cumulant terms and dynamic fermions 
in the calculation of the Wegner- Wilson loop 
in the stochastic vacuum model 

In this appendix we discuss first some problems arising in the calculation of 
the Wegner- Wilson loop in the SVM (cf. sect. 2.5) if higher cumulants are 
taken into account. Then we outline a possible remedy which may also point 
to a way of including the effects of dynamical fermions in the SVM. We start 
with the replacements (|2.79| ) which allow us to use the cumulant expansion 
( |2.41| ) for calculating W(C). The second cumulant is given in ( |2.80| ): 



K 2 (1,2) = ~F{1,2), (A.l) 
3g 2 

where we set with as defined in ( p. 61] ): 

F(iJ) = F 1414 (XW-X^), (A.2) 




(A.3) 
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With the assumptions 1-3 of the SVM (sect. 2.4), the next nonvanishing 
cumulant is K±, for which we obtain from ( p. 47] ), ( |2.60|) , ( |2.73|) , ( |2.75|) : 



/ 2\ 2 

AT 4 (1,2,3,4) = -2 f^J [^(1,3)^(2,4)9(1,2,3,4) + perm.}. (A.4) 

Here 6(1, 2, 3, 4) = 1 if X« > X® > I< 3 ' > I< 4 ) in the sense of the path- 
ordering on the surface S (cf. Figs. 8, 10) and 0(1,2,3,4) = otherwise. 

We start now again from the expression ( |2.7S| ) for the expectation value 
of the Wegner- Wilson loop and use the cumulant expansion ( |2.41| ) with the 
identifications (|2.79|) . We get then: 



W(C)=exp\jr { -^ [ ... / K n (l,...,n)\, (A.5) 
n\ J Si Js n J 



where 

dXPdXf. (A.6) 



is, Js 

In the SVM as formulated in sect. 2.4 the cumulants for odd n vanish 
(cf. Ass. 3). The lowest contribution in (|A.5| ) is then from n — 2, the next 



from n — 4. Cutting off the infinite sum in (|A.5|) at n = 4, we get with ( A.lj) 
and ( [All) : 



H/(C) 



i^/ si .../ s ^ l( l,2,3,4)} 



4! 

exp 



1 4tt 



2 



21 3 JSi Js 2 

-^vr 4 ^...^(1,3)^(2, 4)6(1, 2, 3, 4) + perm.^A.7) 
Consider now the contribution of the second cumulant in ([A. 7 ) for a large 



Wegner- Wilson loop (Fig. Al with R, T > a): 

I 2 = I I ^(1,2). (A.8) 
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For fixed integration point (1) on S the integration over (2) will give signifi- 
cant contributions only for a region of radius ~ a around (1) since the func- 
tions D(-Z 2 ) and D^-Z 2 ), where Z = (X^-X (2 ^), are assumed to fall off 
rapidly with increasing Z 2 . From ( |A.3| ) we see that the Dx-term contributes 



as a total divergence in .F(l, 2). Thus we can apply Gauss' law in 2 dimen- 
sions for it to transform it to an integral over the boundary dS = C. In this 
way we find a contribution from the Di term of order a A /R 2 , a 4 /RT, a A /T 2 . 
From the D-term in ( A.3|) we get a factor oc a 2 from the integration over (2) 



in ( |A.8|) . Then the integration over (1) is unconstrained and gives a factor 
proportional to the whole area of S: 

I 2 ocRTa 2 + o(a\a^,a^). (A.9) 
V R T J 



Putting in all factors from ( |A.3|) and using the explicit form (|2.65|) for the 



function D(—Z 2 ) gives for T — > oo and R ^> a the result: 

1 4vr 2 



2! 3 



I 2 = RTa (A. 10) 



with a as given in ( 2.83|) 



We turn next to the contribution of the 4th cumulant in ( A.7| ) and study 
the integral 

h= f ... f .F(l, 3)^(2, 4)0(1, 2, 3, 4). (A.ll) 

JSi Js 4 

Consider again R,T 3> a and fix the integration point (1) (Fig. A2). Then 
the short-range correlation of the field strengths requires (3) and (1) as well 
as (4) and (2) to be near to each other, i.e. at a distance ~ a. The function 
0(1, 2, 3, 4) requires (1) > (2) > (3) > (4) in the sense of the path ordering on 
the surface 5*. Using the fan-type net as indicated in Fig. 8 for the application 
of the non-abelian Stokes theorem with straight line handles from the points 
X^> in the surface to the reference point Y we see the following. The path- 
ordering function 0(1,2,3,4) restricts (2) to be in the hatched sector of S 
in Fig. A2. In general this does not restrict (2) to a region close to (1). On 
the contrary, (2) can vary freely over a triangle of area ~ a.L, where L is 
of order R or T or in between. Thus, the integration over (2) in ( A.ll|) will 



give at least a factor oc (aT) if the hatched triangle has its long sides in X$ 
direction (cf. Fig. A2). The integrations over (3) and (4) in ([A.ll ) should 



61 



give factors of a 2 each. The integration over (1) will give a factor RT. Thus 
we estimate 

I 4 oc RT.a 4 .aT oc T 2 . (A. 12) 

This is very unpleasant. The 4th cumulant gives a contribution which domi- 
nates over the one from the second cumulant for T — ► 00. There is no finite 



limit for T — > 00 in the expression ( [2.771) . The quark-antiquark potential 
comes out infinite. 

This problem was recognized in ref. |52| and eliminated by hand making 



an additional assumption: that all but the leading powers in the quotient of 
the correlation length to the extension parameters R, T of the loop could be 
neglected. Another simple cure of the problem would be to postulate instead 
of Ass. 3 (cf. (2.73)-(2.75)) a behaviour of the higher point correlation 
functions of the shifted gluon field strengths which gives precisely zero for 
the cumulants K n with n > 2 in (A. 5). In our opinion this would be an 
unsatisfactory solution, since the model would then only work for a particular 
fine-tuned set of correlation functions whereas generically the above problem 
would remain. 

We think that the origin of these difficulties in the SVM is the assumption 
1, which states that the correlator should be independent of the reference 
point Y for arbitrary Y. But why should the field strength correlation ( 2.6Up 



(cf. Fig. 5) be the same if we use a straight line path Cx + Cx> to connect 
X and X' or a path which runs on a loop behind the moon? We will replace 
Ass. 1 by a milder one: 



Ass. 1': in (|2.60| ) is independent of the reference point Y and of 



the curves Cx and Cx> if the reference point Y is close to X and X': 

\X-Y\ &a', \X'-Y\ &a', 
where a' is of order of a. 



Now we try to calculate the expectation value of the Wegner- Wilson loop 
fl2.76p using only this weaker hypothesis. We start from the rectangle S of 
area RT and insert a smaller rectangular loop C\ with sides R — 2a', T — 2a'. 
On C\ we choose N\ reference points Y"i, and we partition the area 
between C and C% in Ni plaquettes Pi, ...,P/Vi of size ~ a' 2 (Fig. A3). We 
can now apply the non-abelian Stokes theorem. We start from Y Q on C and 
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construct a path equivalent to C in the following way: From Yq to Y±, then 
in a fan-type net over the plaquette Pi with Y\ as reference point, etc., until 
we arrive at Y^ 1 = Y\ from where the path runs back to Y Q . According to 



( |2.36|) we get then: 



W(C) = l -TT{V^ Nl V{P Nl )V Nl ^ l .V 2tl V{P 1 )V lfi ), (A.13) 



where 



V(P 5 ) =P expi-ig f G 1 ^Y j ,X^;C x(J) )] (A.14) 

jPj 

and Vij are the connectors from Yj to Yj on straight lines. 

• Ass. 4: Now we will make a mean field-type approximation and re- 
place the path-ordered integrals over the plaquettes Pj by the corre- 
sponding vacuum expectation values: 

V{P S ) - <y{P s )) ■ 11. (A.15) 

This is similar in spirit to the "block spin" transformations considered 
in [j55|] . For the r.h.s. of (|A.15| ) we can use the cumulant expansion and 
assumptions 1', 2 and 3 of the SVM. Here the reference point Yj is never too 
far away from We get with H the unit matrix in colour space 

(V(P 1 )) = l-ex P {f:tl9r f K n (l,...,n)} 

71=1 JP J.l J "j,n 

= 11 ■ exp { - \Pj\a[l + 0{^a ■ a ■ a')}}. (A.16) 

Here we cut off the cumulant expansion at n = 4 and use the results and 
estimates (|A.7|) - (|A.12|) , but now for each plaquette Pj instead of the whole 



surface S. We see that the correction terms from the 4th cumulant are now 
of manageable size since (cf. Q2.70| ), ( |2.85|) ) 



aaa' = aa 2 = 0.56 (^j = 0(1). (A.17) 

Then (hopefully) the factorials 1/n! in the cumulant expansion will lead to 
small corrections from higher cumulants. 
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We can now insert ( |A.15| ), ( |A.16| ) in ( |A.13| ) and get 



W(C) ~ -Ti(V , Nl (y(PN 1 ))V NuNl . 1 ...V 2A (V(P 1 )}V lfi } 

= \Tr(V 0;Nl V Nl)Nl ^...V 2>1 V 1>0 ) -exp{ - ]T (P^l + O^)]} 

d 3=1 4 - 

= ^(^^-1.-^1) exp { - - \Si\)a[l + O(^)]}. (A.18) 



Here we used also the cyclicity of the trace and Vi oV^tvi = U (cf. ( |2.24|) ). The 
product of the remaining connectors in ( |A.18j ) gives just the Wegner- Wilson 
loop of the curve C\ = dSi in Fig. A. 3 and we obtain: 



W(C) = W(d) exp { -(\S\ - \Si\)a 



1 + 



aaa 



4! 



(A.19) 



The procedure can easily be repeated by inscribing a rectangle 5*2 with bound- 
ary C*2 = dS% in Si etc. The final result of this iterative procedure obviously 
is: 

' ' ' aaa: 



4! 



W(C) =exp j-|S|<7 1 + 

From this we deduce for the "true" string tension in the SVM: 

aaa' 



' true 



a 



1 + 



4! 



a [1 + terms of order 



(A.20) 



(A.21) 



where a is given in ( 2.83|) . Thus we have justified the result for the string 
tension in section 2.5, where we used the second cumulant only. We have 
now relied on the assumption 1' and avoided transportation of field strengths 
to far away reference points Y. The fourth cumulant is estimated to give 
only a correction of a few percent and the contribution of the 6th, 8th, etc. 
cumulants can be estimated in an analogous way to be even smalleif] Most 
important, we have avoided the unpleasant result ( A. 12 ). 

We will now discuss another question which can be raised in connection 
with the calculation of W(C) in the SVM. Why should we span a minimal 

2 Numerical studies suggest that a' w 3a should be large enough for obtaining the area 
law for the plaquettes Pj in ( A. 16 ) (H. G. Dosch, private communication). We obtain 



then aaa' ~ 1.7 and still correction terms ~ 10 % in (A.21) 
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surface S into the loop C and not use some other, more complicated surface 
S' with the same boundary, dS' = CI From the point of view of applying 
the non-abelian Stokes theorem ( J2.37|) any wiggly surface S' would be as 
good as the flat rectangle S which is the surface of minimal area. However, 
from the point of view of the iterative calculation with the plaquettes, as 
explained in this appendix, an arbitrary surface S 1 is clearly not equivalent 
to S. In (|A.15|) we made the approximation of replacing V(Pj) (j = 1, Aq) 



by its vacuum expectation value. This should be a good approximation if 
the various plaquettes are well separated. Indeed, we would expect that 
then we can perform the functional integral over the variables related to the 
regions in Euclidean space time where the various plaquettes are located in 
a separate and independent way. If, however, two plaquettes overlap or are 
very near to each other, the above approximation will break down. The 
point is now that on the minimal surface S neighbouring plaquettes will 
be at maximal distance from each other. For some arbitrary surface S' with 
wiggles we will inevitably find plaquettes closer together which will make our 
approximation worse. This gives us some rationale to use a minimal surface 
S in the applications of the non-abelian Stokes theorem in the framework of 
the SVM. 

So far our calculations should apply to QCD with dynamical gluons and 
static quarks only. The quark-antiquark potential V(R) in ( |2.84j ) rises lin- 
early for R — > oo. In real life we have, of course, dynamic quarks. If we 
separate a heavy quark-antiquark pair QQ starting from small R we will 
first see a linear potential as in ( p. 84 ) but at some point the heavy quark Q 



(antiquark Q) will pick up a light antiquark q (quark q) from the vacuum 
and form a meson Qq (Qq). The two mesons can escape to infinity, i.e. the 
force between them vanishes, the potential V(R) must go to a constant as 
R — > oo: 

V(R) -> Voo for R -f oo. (A.22) 

Can we understand also this feature of nature in an extension of the SVM? 

Let us go back to the approximation ( A.15p , where we have replaced the 
integral over the plaquette Pj by its vacuum expectation value. More gener- 



ally we can argue that V(Pj) as defined in (|A.14j) is an object transforming 



under a gauge transformation fl2.17|) as follows (cf. ( |2.20| )): 

V(P j ) - U(Y j )V(P j )U~ 1 (Y j ). (A.23) 
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Any approximation we make should respect this gauge property and indeed 
the replacement ( A.15Q does. Now we can ask for a generalisation of (|A.15|) in 



the presence of dynamic light quarks. We have then the quark and antiquark 
variables at the point Yj at our disposal and can construct from them the 
object 

q(Y 3 )q(Y 3 ) (A.24) 



which has the gauge transformation property ( |A.23| ) and is also rotation- 



ally invariant. Let us consider only u and d quarks. Then we suggest as 
generalization of ( |A.15| ) the following replacement: 

Vab(Pj) - {ViP^AB ~ w (Pj)QAJ,a(Yj)qBJ,a(Yj), (A.25) 

where A, B are the colour indices, / = u, d is the flavour index and a the 
Dirac index. Furthermore, (V(Pj))o is as in (A. 15), (A. 16), and w{P 3 ) is a 
coefficient depending on the size of the plaquette P 3 . It can be thought of 
as representing the chance of producing a qq pair from the vacuum gluon 
fields over the plaquette P 3 . (This is inspired by the discussions of particle 
production in the LUND model |54j| ). On dimensional grounds w(Pj) must 
be proportional to a volume, thus we will set 

w(Pj) = c\Pj\ ■ a, (A.26) 

where c is a constant. The idea is that qq "production" should feel the gluon 
fields in a disc of area \P 3 \ and thickness a. 



We can now insert the ansatz ( |A.25| ) in ( |A.13j ). The resulting expression 
for W(C) is of the form: 

W(C) £ ^Tr(V , Nl [(V(P Nl )) 

-w (P Nl )q(Y Nl ) q(Y Nl ) ] V Nl ^-1 

...V2,i[(V(Pi)) - w(P 1 )q(Y 1 )q(Y 1 )] ■ V h0 ). (A.27) 

If we mutliply out these brackets we get terms where we have again the 
Wegner- Wilson loop along C\ (Fig. A. 3) but then also terms where quarks 
and antiquarks at different points Yfc, Y\ are connected. For two neighbouring 
points, for instance, k = I + 1, / this would read: 

...g(r m )[- w (/)g(r m )y m ^(rO]^)- 
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The importance of these terms will increase with increasing w(Pj). Starting 
from ( A.27 ) we can now inscribe plaquettes into the rectangle Si and trans- 



port everything, including the quark variables q(Y k )q(Y k )(k = 1, ...,Ni) to a 
curve C 2 etc. 

To get an orientation we will now make a drastic approximation: For a 
loop with R <C R c , where R c is some critical value to be determined, we 
neglect the dynamical quarks, saying that the w(Pj) factors will still be too 
small. The result for W(C) and V(R) in the range a ~ R <C R c will then be 
as in (pC20|) : 



W{C) =exp(-RTa), 
V(R) = aR, 

(a^R<^R c ), (A.28) 

where we set a truc = a. 

In the other limiting case, R 3> R c , we can expect the qq terms in ( A.27] ) 



to dominate. Indeed, let us start with the loop C with sides R and T of Fig. 
A. 3 and inscribe first plaquettes Pi v ..,P/v 1 of size a' 2 . With the replacement 
( |A.25| ) this gives for W(C) the expression ( |A.27| ). Now we inscribe plaquettes 



of size a' 2 in the loop C\ and parallel transport everything to a smaller loop 
C*2. The new element is that we will now have to parallel-transport also 
the quark variables, but this poses no problem. We will again obtain an 
expression like the one shown in ( |A.27| ) but now for the curve C 2 and with 
increased weight factors w in front of the qq terms. We assume that we 
continue this procedure. Finally, we will obtain an expression like ( |A.27j ) 



but with the terms involving the quark variables at the reference points Yj 
dominating the expression. Thus we set: 

W(C) S ^Ti([-w(N)q(N)q(N)}V N , N ^...[-w(l)q(l)q(l)}V 1 , N ), (A.29) 

where the points 1, N are on a curve C of sides R', T' . This summing up of 
quark contributions from various plaquettes of size a' 2 should be reasonable 
as long as the quarks are inside an area corresponding to their own corre- 
lation length, a x , the chiral correlation length. This can be estimated from 
the behaviour of the non-local qq condensate [86-88] for which one typically 
makes an ansatz of the form (q = u,d): 

(q(X 2 )V 2il q(X 1 )) = (qq)e-^- x ^< (A.30) 
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Here (qq) is the local quark condensate for which we set, neglecting isospin 
breaking (cf. j39| ): 



(m) = \(uu + dd) = -(0.23 GeV) 3 . 



(A.31) 



For the correlation length a x one estimates (cf. fl38|| ): 

a x ^ 1 fm. (A.32) 

It will be reasonable to choose the reference points 1, ...,N on C also such 
that the distance between neighbouring points is a x . Then 

R! 



r 

N 



R-2a 
T -2a 



2V + 2R' 2{T + R) 



(A.33) 



X X 

We will now estimate W(C) from flA.29 ) as product of the expectation 
values of the nonlocal qq condensate. We set 



w(l) = w(2)... = w(N) = c ■ aa 2 x , 



(A.34) 



where c should be a numerical constant of order 1. Furthermore we have 
from ( |Oq ) 

(QAJ,a{j)( V h3-l)AB QBJ',a'{j ~ 1)) 



exp 



\Xj — XjA 



= —5 fif 8 a , al (qq) 
for \Xj — Xj_i \ = a x . This gives: 



(A.35) 



W(C) 



-.{[-wiN^qf^l^N q fN , aN (N))} 
[-w(l){q h>a2 (2)V 2>1 q fl , ai (l)}}} 
w(l)(qqy A 



4e 



w(l)(q q) 
4e 



(A.36) 
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For the potential this leads to 
V(R) = 



1 2 4e 
lim -\nW(C) = — ln T , w n 



(A.37) 



Thus we find indeed a constant potential for R — ► oo. 

For intermediate values of i? the potential should change smoothly from 
the linear rise to the constant behaviour. As a crude approximation we will 
assume the potential to be continuous but turn abruptly from one to the 
other behaviour at a critical length R = R c which we define as 



R c = Voo/a. 

Our simple ansatz for V(R) is then as follows (cf. Fig. A. 4): 

aR for R < R r , 



V(R) 



aR r for R> R r . 



(A.38) 



(A.39) 



How does the numerics work out taking ( |A.37| ) and ( |A.38| ) seriously? 
We can estimate the constant V<x> crudely from the mass difference of two 
.B-mesons and the T resonance 



Voo = 2m B -m r £ 1.1 GeV. 
With a from ( ggg) we find for R c (pQ8|) : 

R r ^ 1,2 fm. 



(A.40) 



(A.41) 

A much more elaborate estimate of the string-breaking distance gives [|89| 



R c < 1,6 to 2, lfm. The lattice calculations of Q suggest an even larger 
value for R c . From (TA~40|) and ( |Q7D , (|02D we get for iw(l) and c: 



"(1) 



4 exp[l — \a x V a 



a x a 



(qq) 

1.2. 



0.42 W 



(A.42) 



Thus the probability factor w(l) comes out as we estimated it on geometrical 
grounds in ( A.26[) . 

With these remarks we close this appendix. Of course, much more work 
is needed in order to decide if our simple ansatz for incorporating dynmical 
fermions in the stochastic vacuum model is viable or not. 
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Appendix B: The scattering of gluons 



In this appendix we will discuss the contribution of gluons to the scattering 
amplitude of a general parton reaction, generalizing fl3.50|) : 



G(l)+G(2) + ... + q(k) + q(k') + ... -> G(3)+G(4) + ... + q(l) + q(l') + .... (B.l) 



Here, as in all of sect. 3, we use the convention that partons with odd (even) 
number are moving fast in positive (negative) x 3 -direction. Consider thus 
the transition of a right-moving gluon: 



where a^3 (1 < 01,3 < 8) are the colour indices, P±, P3 are the momenta with 
Pi + ,P% + — > 00, and €1.3 are the polarization vectors which satisfy: 



In the high energy limit the vectors e 13 are transverse with corrections of 
order |pri,3|/-P+i,3 which we can neglect. We will argue that such a gluon 
in a high energy soft reaction is equivalent to a quark-antiquark pair with 
the same quantum numbers in the limit of the q and q being so close to 
each other in position space that their separation cannot be resolved in the 
collision. From the contribution of such a qq pair to the scattering amplitude 
we will get in the above limit the gluon contribution. 

We start by constructing wave packets similar to the mesonic wave packets 



G(l) = G(P U e 1 , ai ) - G(3) = G(P 3 , e 3 , a 3 ), 



(B.2) 



Pl6i = 0, 

P 3 e 3 = 0, 
|e x | = |e 3 | = 1. 



(B.3) 



of (gg§) 




(B.4) 



where pj is as in ( [3.67D , (|3.68|) and 




(B.5) 
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We require the functions h? to satisfy: 

(tf(C,|Pr|))* = ^'(C,|Pt|), 

^(C,Ipt|) = ^'(i-C,Ipt|), 

J d 2 p T | o 1 dC2C(l-C)l^'(C,|PT|)| 2 = l. (B.6) 
This gives us the normalization of the states (|B.4|) as 

= <': ■ c,6,„.„ (2^) :: 2/^ :! (l^ - !>,•). (B.7) 



The qq states (|B.4|) have the same transformation properties under a 
parity (P) charge conjugation (C) and time reversal (T) transformation as a 
gluon state: 

U(P)\qq; P°, P, e, a) = -\qq; P°, -P, e, a), (B.8) 
U(C)\qq; P°, P, e, a) = -\qq; P°, P, e, b) • ^Tr(A b A^), (B.9) 
y(T)|gg; P°, P, e, a) = -\qq- P°, -P, e*, b)~Tr(\ b \*). (B.10) 

Here t/(P), i7(C) are the unitary operators, V"(T) is the antiunitary operator 
representing the P, C and T transformations, respectively, in the Fock space 
of parton states. 

As in ( |3.71| ) we define the wave functions in transverse position space and 
longitudinal momentum fraction: 

^'(CYt) := V^Ca-C)^ / d 2 p T exp(tp T ■ y T )h j ((, |p T |). (B.ll) 

Here we also define for the right (left) movers the wave functions in y-,yT 
space (y + , yr space) as: 

^(y-,y T ) ■= -7=^ f^C ^(C,yr)exp[-lp i+ (c- \)y-\, 



'Air 

U odd), (B.12) 



/4vr 

U even )' (B.13) 
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The normalization condition ( |B.6| ) implies: 



12 -1, 



dV /'dCI^'^y^P 

JO 

J dy„J d 2 y T \ft(y„,y T )\ 2 = 1 (j odd). 

J dy + J d 2 y T \^(y + , y T )\ 2 = 1 (j even). (B.14) 

To realize the condition that the qq pair acts like a gluon, we require 
for right (left) movers that they have similar longitudinal momenta and that 

their wave function in the relative q-q coordinates, y_, yx(y+, yr) is nearly a 
5 function. To be concrete, we require: 

^(CYt)^O only for |C-^|<6>, (B.15) 

where 

< fo « \, (B.16) 

and for j odd: 

<^(2/-»yr)^0 only if \y_\ ~ — — and |y T | < a. (B.17) 



For right movers we replace plus by minus signs in (|B.17f ). Any qq wave 
packets with these properties should then look identical to a gluon for an 
observer in the femto universe with regard to "soft" scatterings. 

Now we replace G(l) and (2(3) in reaction ( |B.1| ) by the qq wave packets 
( |B.4| ). According to the rules derived in sect. 3.5 the scattering of the qq 
system 

g(l)g(l') - g(3)g(3') (B.18) 
gives the following factor in the S'-matrix (cf. ( |3.58| ) (ft.59|) ): 

which still has to be integrated over the wave functions ( |B.4[ ). In this way 
we obtain for the gg-pair: 
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5^(3,1) = J d 2 p' T J\cj d VjV 



^ 3 (C',IptI)^(C,Ipt| 



<S g+ (3, l)5 ff+ (3',l'), 



(B.19) 



5^(3, 1) = / dz_ | dy_ | d 2 2T | d 2 y T JP^P^e* 3 ■ e x 
exp[^(P 3+ - Pi+)z- - i(P 3 - Pi)tZt] 

( ^ 3 *(z/-,yT)^ 1 (y-,yT) 

1 11 

-Tr[A a3 F_(oo, 2_ + z r + -y r ) 
A ai 1/i(oo, z_ - z T - ^y r )]. 



(B.20) 



With our assumptions (B.17| ) Lp 3 *(y_,y T )ip 1 (y_,y T ) acts like a 5-function at 
y_ = 0, = and we get: 



1) — yP 3+ P 1+ €3 • e x 

J dy_ /" d 2 y T <^ 3 *(y_,y T )^ 1 (?/_,y T ) 

| dz_ | d 2 z T exp[^(P 3+ - Pi+)z- - *(P 3 - Pi)t • Zr] 



-Tr[A a 3y_(oo, z_, z T )A ai Vi(oo, z_, z T ) 



(B.21) 



An easy exercise shows that 



^Tr[A 03 V_(oo,z_,z r )A ai Vi(cx),z_,z T )] = V_(oo, z_, z T ) 



03,01 j 



(B.22) 
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where V_ is the connector analogous to ( |3.35| ) but for the adjoint represen- 
tation (cf. (ggp ): 

i f°° 

V-(oo,z_,z T ) = P{exp [- -g / dz + G a _{z + , z_,z T )T a ]}, 

Z J — oo 

{T a )bc = -fabc- (B.23) 

% 

The f a bc are the structure constants of SU(3) C . Inserting ( |B.22|) in (|B.21| ) 
gives 

S qq (3,l) =Jdy-J dV ^ 3 *( 2 /_,y T )^ 1 ( Z /_,y r )e^ eijl 5 G+ (3,l), (B.24) 
where 

S G+ (3, 1) = V^P^ / Ck_ / d 2 Z T 

exp [^(P 3+ - Pi+K - *(Psr - Pit) • z r ] 

V_(oo,2_,z T ) a3jai . (B.25) 

In a general scattering reaction ( |B.1| ) the factor S qq (3, 1) ( |B.24| ) has to 
be inserted with other factors S q ,S q ,... and then integrated over all gluon 
potentials as explained in sect. 3.5. We note that S qq (3, 1) in ( |B.24p factorizes 
into Sq + (3, 1) ( |B.25|) times the overlap of the internal wave functions of the 
incoming and outgoing qq pairs: 

J dy-J dV ^*{y^j T )(p\y^y T ). (B.26) 

This means that the qq pair will come out with some distribution in total 
momentum P3 and polarization vector €3 = e\ but always with internal wave 
function ip 1 (y_,yx). With the conditions (|B.15|) -( |B.i7[ ) <^ 1 (?/_,yr) leads to 
a "permissible" internal wave function for a qq pair of momentum P3, to be 
regarded as a gluon of momentum P3 by our observer in the femto universe. 
Indeed we have from ( p,12|) : 



v\y-,YT) = ^j\ccp\c,y T )eM-lPi+(C-~)y-) 

I P, , rl/2 1 j 



An J$ 



d£' + y T ) exph^Pa+C • y_](B.27) 
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Here we define the internal wave function of the outgoing qq pair as 
1 

and 




^ (I + Yt ) = J^rV ( ^ + ^C, Yt ] (B.28) 



e ± = ±^ • J- (B.29) 

The condition ( |B.15[ ) for if 1 guarantees a similar condition for ip' 3 if P 3 _|_/P 1+ 
is of order 1, as we will always assume. Thus the integration limits £± in 
( |B.27| ) can be replaced by ±~ and also the normalization conditions ( |B.14| ) 
can easily be checked for ip' 3 . 

To summarize: In this appendix we have shown that suitable qq pairs 
which are indistinguishable from gluons for our observer in the femto uni- 
verse scatter as entities in a soft reaction. Their internal wave function in 
momentum space is modified, but in a way not observable in a soft reaction. 
Thus we can consider <S G+ (3, 1) in (|B.25|) as the scattering contribution of 



a right-moving gluon in the transition ( |B.2| ). We quoted this result already 
in (|3.64j ) in sect. 3.5. For left moving gluons we just have to exchange 
everywhere + and — components. 



Appendix C: The scattering of baryons 

In this appendix we discuss high energy soft reactions, in particular elastic 
reactions involving baryons and antibaryons. We represent baryons by qqq 
wave packets: 

\Bj{Pj)) = -^^y J dfi U{f\s\C\\tr) ■ 

ca^aMp], f\ s\ A 1 ), q(pl f, s 2 , A 2 ), q(p 3 , f, s 3 , A 3 )). (C.l) 

Here /*, s l , A 1 (i = 1, 2, 3) are the flavour, spin, and colour indices of the 
quarks and €abc ( e i23 = 1) is the totally antisymmetric tensor. For the 
momenta we set for right-moving baryons (j odd, Pj + — > oo): 

I', - = CI).. 

(* = 1,2,3) (C.2) 
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and for left-moving baryons (j even, Pj_ — > oo): 

v)- = CP,-, 



3 

i = 1,2,3). (C.3) 



The integral with the measure d/z stands for 

3 3 



j j /ndv T 5 2 (Ep^)- 

i=l i=l 
rl rl rl ^ 3 

/ / / nw-Ec). (c.4) 



The flavour and spin of the baryon states (|C.1|) is, of course, determined 
by the functions h? which must be totally symmetric under simultaneous 
exchange of the arguments (/*, s l , ( l , p l T ) for i = 1,2, 3. In the following we 
will collectively set (/, s, (, p-r) = ol and q(a, A) = q(j)j, f,s,A). We have 
then 

h j (a,(3,-f) = h j ((3,a-f) = h j (a } -f,(3). (C.5) 
\B {P 3 )) = -^-^ / d^ h^a,P, 7 )e ABC \q(a,A)q(P,B)q( 1 ,C)).(C6) 

The normalization condition 

(BiWmPi)) = (2vr) 3 2 J R ? ° ( 5 3 (P;. - P,) (C.7) 
requires h? to satisfy: 

d// 4C 1 C 2 C 3 ^'*(«,/?,7)^'(«,A7) = 1 

no summation overj). (C.8) 



We define the wave functions tp J in transverse position and longitudinal mo- 

% f o r B i 



mentum space and the transition profile functions wE for Bj — > Bk as: 
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^(f\ s\ 0, 4) := tA^ / f[ d "PT ■ Pt) " 

v / i=l i=l 

6 • (CCCT 2 exp(* £ p*r ■ 4M«\ a 2 , a 3 ), (C.9) 



-I r l rl 3 



<K,x 2 ,x 3 ):= / / / IldC^l-EO 

J0 J0 i=i i=i 

E <^*(A C 1 , x^; / 2 , a 2 , C 2 , x 2 ^; f, s 3 , ( 3 , x 3 ) 

/(/\ a 1 , C 1 , x^; / 2 , , 2 , C 2 , x 2 .; / 3 , S 3 , ( 3 , x 3 ). (CIO) 

The symmetry relations for h ( |C.5| ) and the normalization condition ( |C.8| ) 
imply: 

w k,j( x Ti x Ti' x T) = w k,j( x T^ :K Ti x T) = W k,j( :x -Ti x r, X T ), (C.ll) 



/ J| d 2 4 5 2 (x^ + 4 + 4)™jj(4> 4, x 3 .) = 1 

"* i=l 

(no summation over j). (C.12) 
As a concrete scattering reaction let us consider meson-baryon scattering: 
Si (Pi) + M 2 (P 2 ) -> P 3 (P 3 ) + M 4 (P 4 ). (C.13) 
From the rules of sect. 3.5 we get 



S fl = (B 3 (P 3 ),M A (P A )\S\B 1 (P l )M 2 (P 2 )) = 

dfi' h 3 *(a',(3', i)t A 'B>c I d/i h l {a, (3^)e AB c 



6(2tt) ( 

(S q+ (a, a) A > A S q+ (f3\ P)B'EtS q +{i, j)cc 

5 A > 2A J A > 4A4 S q „(4, 2) AaA2 S^(4', 2')a^) g - (C.14) 
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After some straightforward algebra we get: 
S fi = 5 fi + i(2ir) 4 5(P 3 + P±-P 1 -P 2 )T fi , 

3 

i=i 



where 



n c*o 



i=l 



1/2 



^ / d^V/ dV jf'dC ^(C',Pt)^(C,Pt) 

2(P2-) 2 [C'(i-OC(i-C)] 1/2 

3 3 • y d 2 6 T exp(iq T • br) / n(da; l _d 2 4)<5 (l>-) ^ (j^) 



dy+d 2 y T exp {« ^ 



i=i 



-p 1+ (C ,1 -CR-(pt-Pt)-x^ 



({V-ioo^i, ^b T + x^) AM \/_(oo,x 2 , ^b T + x^) B / B 

1 i 

V_(oo,a;i, -b T + x^)cc g^A'B'c £abc 

^Tr[V + (y + , oo, -^b T + iy T )Vj (0, oo, -^b T - ±y T )] - 1}) G , 

(C.15) 

(C.16) 



— (Pi — Ps)t- 
Now we make the transformation of variables 

2 . 



x 



y+ 



Pi+ 

2 



-x 



-y+ 



(C.17) 



and use P i+ — > oo, P 2 _ — > oo. With the same arguments which led us from 
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(157770 to ( 157750 we get 



7^ = -2zs /" d 2 6 T exp(iq T • b T ) / JJ d 2 x^ <5 2 (x^ + x| + x 

^i( x t ; x t> x t) / d 2 y T w™ 2 {y T ) 

({VI (oo, 0, ib T + x^) j4 / A Vl(oo, 0, ^b T + x 2 .)^ 



V-(oo, 0, -b T + x|.) c ./ C -e yl '_B'c'eABC 
2 o 

-Tr^O, oo, -^b T + iy r )^(0, oo, -^b T - iy T ] - 1}> G . 

(C.18) 

Here w^ivr) is the transition profile function for the mesons as defined in 
( |3.72| ) and is the corresponding function for the baryons (cf. ( |U.1U| )). 

In the next step we follow |p9| , |f45f and use relations which are valid for 
any 3x3 matrix: H = (Hab)'- 

Ha'A"Hb'b"Hcc" ' €a"B"C" — detH ■ €a>b'c, (C.19) 



detif • eA'B'c'tABC — 

Ha'aHb'bH C 'c + Ha'bHb'cHc'a + H A 'cHb'aHcb 

—Ha'bHb'aHcc — Ha'aHb'cHcb — Ha'cHb'bHca- (C.20) 

We take H equal to the antiquark line integral at the central point of the 
baryon in transverse space: 

# = VL*(oo,0,^b T ). (C.21) 
As a Si7(3)-connector V* satisfies 

detV!(oo,0,-b T ) = I. (C.22) 



This is easy to prove. From ( |C.19| ) we see that det V* is the connector in 



the singlet part of the product of three SU(3) antiquark representations: 
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3x3x3. But for the singlet representation the connector equals 1 since we 
have to set T a = in (ggg ). 



In the following we will use as shorthand notation 

^(i)=V_(oo,0,ibr + 4), 
(i = 1,2,3), 

y(0)=V_(oa,0,~br). (C.23) 

With this we get for the ggg-contribution to the integrand in the functional 
integral in (|CJj) using (p^Ol) to (|02l ): 



W|(-bT,4,x^x|,) : = -y(l) A 'AV(2) B / B y(3)cc e^ceABc 
2 o 

= ~ (Tr[01)0(O)] • Tr[02)0(O)l • Tr[03)0(O)l 
6 L 

+ ^[01)0(0)02)0(0)03)0(0)] 
+ Tr[01)0 (0)03)0(0)02)0(0)] 

- Tr[01)0(O)02)0(O)] ■ Tr[03)0(O)] 

- Tr[02)0(O)03)0(O)] ■ Tr[01)0(O)] 
-Tr[01)0(O)03)0(O)] ■ Tr[02)0(O)]} . (C.24) 

As for the meson case (cf. (|3.82|) ) we will now add suitable connectors at ±oo. 
In this way all the traces in ( |C.24j) become closed light-like Wegner- Wilson 
loops. To give an example: the term 

Tr[01)0(O)02)0(O)] 

should then be read as the loop in the hyperplane x_ = in the limit T — > oo 
which connects the following points (x + ,x_,x r ): 

(r,o,ib r ), 

(-T,0,ibr), 
(_ T)0 ,ib r + x^), 
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(T,0,-br + 4), 

(-T,0,ib r ), 
(_ T)0 ,ib T + x^), 
(T,0,^b T + x^), 

(T,0,~b T ) (C.25) 

on straight lines in the order indicated. 

Inserting (|C.24|) in ( |C . 1 8|) and denoting the mesonic Wegner- Wilson loop 
as defined in ( |3.83| ) by W* 1 , we get finally for the T-matrix element of baryon- 
meson scattering: 

Tfi = —2is J d 2 bx exp(zq^ • by) 
f 3 

/ Yl d 2 Xy S 2 (x}p + x^ + x|,)w^ 1 (xy, 4> x|) 
/ d 2 ?/T ^(yr) 

(W* Qb r ,4,x3. l x§.) W M (-^br.yr) - . (C.26) 

This formula is the starting point for the evaluation of the baryon-meson 
elastic scattering amplitude: One can now apply the Minkowskian version 
of the SVM to calculate the functional integral ( }g in an appropriate way. 
Then one has to fold the result with the profile functions of the mesonic and 
baryonic transitions. At the present state one has to make a suitable ansatz 
for these profile functions. 

For the case of right-moving antibaryons in an elastic reaction we just have 
to substitute the loop factor by which is obtained by replacing the 
quark connectors V_ by the antiquark connectors V* and vice versa in ( |C.23| ) , 
( |C.24j ). In an equivalent way we can get from by reversing the 
directional arrows on all Wegner- Wilson loops obtained in the way discussed 
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above from ( |C24j ). For left-moving baryons and antibaryons we have to 
exchange + and — components. 

For the further treatment of scattering amplitudes involving mesons, 
baryons and antibaryons, for many results and a comparison with experi- 
ments we refer to |59fl) [EH], |62 
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Figure Captions 



Fig. 1: The lowest order diagram for the Drell-Yan reaction ( |1.9|) in the 
QCD improved parton model. 

Fig. 2: The schematic behaviour of the vacuum energy density e(B) as func- 
tion of a constant chromomagnetic field B according to Savvidy's cal- 
culation (eq. (|2.2|)). 

Fig. 3: A "snapshot" of the QCD vacuum showing a domain structure of 
spontaneously created chromomagnetic fields. 

Fig. 4: The QCD- vacuum according to Ambj0rn and Olesen |33| (a). The 
ether according to Maxwell |38| (b). 

Fig. 5: Points X, X', Y in Euclidean space time and curves Cxi Cx 1 running 
from X to Y and X' to Y, respectively. 

Fig. 6: Curves in Euclidean space time: C\ going from X\ to X 2 , C 2 from 
X 2 to X 3 , C from X to Y, and G from Y to X. 

Fig. 7: A surface S with boundary C = dS. A plaquette with corner points 
Zi..., Z 4 and boundary formed by lines u = const, v = const.. 

Fig. 8: A reference point Y on S and the fan-type net with center Y spanned 
over S. 

Fig. 9: The correlator function 47r 2 KG2-D(— Z 2 ) (cf. fl2.61J) ) as function of 
\Z\. The dashed line is the lattice result of [|0J] with the arrows indi- 
cating the range where these results are considered reliable. The solid 
line corresponds to the ansatz ( |2.65| ) with a = 0.35 fm (cf. Fig. 5.1 of 

Fig. 10: Rectangular Wegner- Wilson loop in Euclidean space-time in the 
X\ — X4 plane. The linear extensions are R in X\ direction and T in 
X4 direction. Y is the reference point used in the application of the 
non-abelian Stokes theorem. 

Fig. 11: The t-channel (a) and w-channel (b) exchange topologies for the 
diagrams describing quark-quark scattering. 
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Fig. 12: Projection of the world lines of the quarks 1(2) moving at high 
velocity in positive (negative) x 3 direction onto the x° — x 3 plane in 
Minkowski space. The non-abelian phase factors V_ in fl3.35| ) (V + in 
( |3.37|) ) are the connectors taken along lines x_,x T = const. (x + ,x T = 



const.). 

Fig. 13: Two light-like lines on which the associated string operators V± in 
(|3.45|) are evaluated. Their correlation function governs quark-quark 
scattering at high energies. 

Fig. 14: Sketch of reaction (|3.50|) in the overall cm. system 

Fig. 15: The light-like Wegner- Wilson loop in Minkowski space-time, C + 
consisting of two light-like lines in the hyperplane X- = and connect- 
ing pieces at infinity. In transverse space the centre of the loop is at 
Yt, the vector from the antiquark to the quark line is zy (from A to 
B). 

Fig. 16: The projections of the two lightlike Wegner- Wilson loops C + ,C_ 
occurring in the definition of W + (\h T , x T ), W-(— |b T , yr) (cf. Q3.871) , 
(|3.83|) ) into transverse space. The points marked q+,q+ (g_,g_) corre- 
spond to the projections of the quark and antiquark lines of C + (C_). 

Fig. 17: The curves C + and C_ along which the path integrals W± in 



(|3.100|) are taken. The mantle of the pyramid with apex at the ori- 
gin of the coordinate system and boundary C + (C_) is V + (V-), the 
basis surface S + 

Fig. 18: The relation between the total cross section a tot and the slope 
parameter b for proton-proton and proton-antiproton scattering. The 
dotted line is the prediction from Regge theory. The prediction of the 
calculation for soft high energy scattering in the stochastic vacuum 
model is that the data points should lie in the area between the full 
lines. In essence this is given by (|3.108|) , (|3.109|) with an uncertainty 



estimate from different assumptions for the proton wave functions (cf. 
151). 



Fig. 19: A proton-proton collision at high energies in the parton picture. 
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Fig. 20: A quark and antiquark traversing a region of chromomagnetic field 
(a). An electron and a positron in a storage ring (b). In both cases we 
expect the emission of synchrotron radiation to occur. 

Fig. 21: Sketch of a "colour domain" in Minkowski space and of the world- 
line of a quark from a fast hadron moving through it. 

Fig. 22: A quark moving in 3-direction in a transverse (in 1-direction) chro- 
momagnetic field of strength gB c and picking up a transverse momen- 
tum (in 2-direction) of magnitude px over a length l e g. Here px is the 
mean transverse momentum of quarks in the hadron. 

Fig. 23: The |k^| distribution for direct photons emitted at cm. rapid- 
ity y = in p — Be collisions at 450 GeV incident proton momentum 
from ||72|| . The normalization is according to a private communica- 
tion by H. J. Specht. The background of decay photons is subtracted. 
The dash-dotted line gives the expected yield of photons from hadronic 
bremsstrahlung, the dashed lines show the upper and lower limits in- 
cluding the systematic errors in the shape of the decay background 
and the bremsstrahlung calculation (cf. |72"| ). The lower (upper) solid 
line is the result of the calculation for synchrotron photons ( (|4.2| ff.) 
with Z e ff = 20 fm (Z e fj = 40 fm) added to the spectrum of hadronic 
bremsstrahlung (cf. |66fl). 

Fig. 24: A nucleon interacting by emission of a quasi-real photon. 

Fig. 25: The data for the electric form factor of the neutron G^(Q 2 ) from 
refs. [7E[| . Dash-dotted line: our naive "synchrotron" prediction 
oc (<3 2 ) 1/6 (|4T3|) normalized to the data at Q 2 = 5 far 2 . Dashed 
line: the slope of G^(Q 2 ) at Q 2 = as deduced from thermal neutron- 
electron scattering |77| . Full line: the ansatz ([4.15|) . 

Fig. 26: The ratio G P E /Gx> of the electric fom factor of the proton to the 
dipole fit versus Q 2 . The data points are from various experiments as 
summarized in jPJ. The solid line corresponds to the ansatz ( |4.17| ), 

Fig. 27: Annihilation of a qq pair with production of a virtual photon 7* 
in a colour domain. Here q and q come from two different hadrons hi 
and h 2 , respectively. 
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Fig. Al: The relevant integration region for the integral ( |A.8| ). The point 



(1) runs freely over the surface S, the point (2) is constrained to a 
distance ~ a from (1). 



Fig. A2: The relevant integration region for the integral ( |A.11|) . The point 



(1) runs freely over S. The points (1), (2), (3), (4) are ordered in 
the angle as seen from the reference point Y due to the path ordering 
function 9(1,2,3,4). The points (1), (3) and (2), (4) must be at a 
distance ~ a to each other. 

Fig. A3: The rectangle S, C = dS and N\ reference points Yi, Y Nl on the 
curve C\. The area between C\ and C is partitioned in N\ plaquettes 
Pi, ■■, Pni- 



Fig. A4: The potential V(R) as defined in (|A.39|) with R c the "string break- 
ing" radius. 
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